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Preface 

This report is concerned with both the theory and details associated 
with the evaluation of large reflector antennas. The work reported herein 
was motivated by this Laboratory's need to both understand and improve 
the performance of the 16-foot radio telescope at the Millimeter Wave 
Observatory, Mt. Locke. Indeed, the heart of this report, the chapters 
on astigmatism, gain measurements, and the case study of the 16-foot 
antenna, describe the fundamental properties of the 16-foot antenna, both 
before and after improvements. 

The earlier Technical Report on the evaluation of the 16-foot antenna, 
TR No. NGL-006-69- 1 , "Calibration Program for the 16-Foot Antenna, " 
by J. R. Cogdell, was by its own admission preliminary and incomplete. 
The present report, taken as a sequel, presents the complete picture of 
the antenna properties. The whole picture is, we feel, very pleasing: 
the antenna surely ranks with the best in its peer group. 

This report is based upon a dissertation submitted to the Graduate 
School of The University of Texas at Austin in partial fulfillment of Ph. D. 
degree requirements. 

The author wishes to acknowledge the commitment and support of his 
supervising professor. Dr. J. R. Cogdell. Thanks also are due to Mr. 
Archie M. Groom for his able technical assistance with the experiments 
reported herein. 

This work is sponsored by the National Aeronautics and Space Ad- 
ministration under NASA Grant NGL 44-012-006. 
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Abstract 


This report concerns the evaluation of large reflector antennas. 

The major findings include defining the effects of primary reflector 
asitgmatism and its quantitative measurement through diffraction pattern 
measurements. A method is also presented which reveals the frequency 
dependence of antenna efficiency. 

Several basic problems are addressed. One is the focusing of 
reflector antennas and the interpretation of reflector errors from the 
antenna pattern. The second topic is the investigation of the effects of 
primary reflector astigmatism. Evidence is presented that several of 
the largest reflector antennas have astigmatism in their primary re- 
flector. The third topic offers a method with which the antenna effi- 
ciency at one frequency can be inferred from the measured antenna 
efficiency at another frequency. The method given here is much less 
restrictive than the Ruze method and bounds the effects of measurement 
errors. 
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Chapter I 


Introduction 

This paper concerns the evaluation of large reflector antennas. 

The major findings include the effects of primary reflector astigmatism 
and its quantitative measurement. A method is also presented which 
provides new insight into the frequency dependence of the antenna 
efficiency. 

The size and precision of reflector antenna systems have greatly 
increased over the last twenty years. Today, the largest reflector antennas 
are used in radar, space communications, and radio and radar astronomy. 
Large, highly precise antennas are required for these applications because 
the received signals are relatively weak. 

When reflector antennas began to appear in the late nineteen- forties, 
there was no theory available to dictate the required precision of manu- 
facture of such antennas. The criterion which had always been used for 
optical reflectors was that the reflector^" ^ must be within X/8 of the required 
shape. However, different criteria are used to judge antennas than are used 

to judge optical instruments, so a reassessment of the problem was in order. 

[ 2 ] 

John Ruze addressed this problem in his doctoral thesis L . His theory 
is now the accepted standard in the field of reflector antenna evaluation and 
design. Probably the most widely used part of his theory is the frequency 
dependence of the antenna gain which it predicts. 
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The Ruze theory is basically a statistical theory and thus is based 
on certain mathematical assumptions. These assumptions specify the 
behavior of an antenna. by implication. Only recently have antenna engi- 
neers come to realize that the Ruze theory is sometimes inadequate to 
describe all of the performance characteristics of large antennas. 

The Millimeter Wave Sciences Laboratory at The University of 
Texas at Austin operates a 16-foot radio telescope. This instrument 
is used in observations of radio sources, especially the planets at fre- 
quencies up to 140 GHz. It was first located in Austin and worked well 

f3l * • 

at 95 GHz. In the winter of 1966/67 it was installed on Mount Locke 
near Fort Davis, Texas to take advantage of the better meteorological 
conditions there. However, soon after it was installed on Mount Locke 
it was determined that the performance of the antenna fell considerably 
short of what was expected. At 95 GHz no satisfactory focus position 
could be found, and the beamwidths were always broad in at least one 
of the principle planes. In addition, one sidelobe appeared at a level of 
-7 dB when the antenna was adjusted for maximum gain. This paper comes 
largely from efforts the laboratory has made to improve the antenna per- 
formance. 

This dissertation addresses several basic problems. One is the 
interpretation of the reflector errors from the antenna pattern. The 
second major topic is the diagnosis of primary reflector astigmatism 
and its measurement. The third topic concerns the dependence of the 
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antenna efficiency on frequency and offers new insight into this problem 
as well as some useful quantitative results. 

Chapter II is a tutorial discussion of the aperture field method of 
reflector analysis. The important underlying ideas of the method are 
discussed with an emphasis on geometric optics. The Fourier trans- 
form relations between the aperture fields and the far zone fields are 
then presented. This chapter may be skipped without harming the con- 
tinuity of the presentation. 

Chapter III is a critical review of the current literature on antenna 
tolerance theory. Particular emphasis is given to the Ruze statistical 
theory, which is widely used by antenna engineers today. In addition, 
the interpretation of antenna patterns is discussed. 

Chapter IV is the presentation of a reflector astigmatism model. 
Several model properties are presented which relate to the far zone 
radiation pattern. These properties are useful in determining whether 
astigmatism is present in a given antenna and also in measuring it 
quantitatively. Evidence is presented through one of these properties 
that astigmatism appears in at least two other large reflector antennas. 

Chapter V gives a method of inferring antenna efficiency at one 
frequency from a measured efficiency at another frequency. If the 
measurement frequency is the higher of the two, the error bars on the 
inferred efficiency are usually less than the measurement error bars. 
Thus, a single antenna efficiency measurement usually serves to define 
the efficiency at all lower frequencies. This method is then extended so 
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that two antenna efficiency measurements can be used to infer the 
efficiency at the third frequency. The proofs for these two methods 
are given in Appendix B. 

Chapter VI is a report of the latest work done in the calibration 
program for The University of Texas 16-foot antenna. It provides a 
concrete example of the use of the astigmatic phase error model in 
Chapter III. A. careful efficiency measurement at 134 GHz is reported. 
The efficiency of the antenna at other frequencies of interest is then 
inferred from a method given in Chapter V. 



Chapter II 


Analysis of Reflector Antennas 

A. Introduction 

At microwave frequencies and above, wavelengths become so short 
that conventional array antennas become impractical. However, short 
wavelength is a desirable feature for reflector antennas, so reflectors 
are more widely used at shorter wavelengths. A reflector antenna system, 
which consists of a small antenna and a reflector, is a simple device that 
focuses electromagnetic radiation in a particular direction. The purpose 
of the small antenna, called a feed, is to launch a wave in the direction of 
the reflector, as shown in Figure 1. The feed usually exhibits a broad 
radiation pattern. At the reflector the shape and direction of the waves 
from the feed are transformed. In particular, they are focused so that the 
radiation of the system is much more concentrated. 

The following section states the general philosophy behind the aperture 
field method of reflector analysis. This is followed by a discussion of geo- 
metric optics. In the final section the Fourier transform relations between 
the far zone fields and the aperture fields are presented. 

B. The Aperture Field Method 

One might analyze a reflector antenna as a boundary value problem. 

However, this is a general mathematical method that yields little insight 

into the actual operation of the antenna. Engineers have traditionally used 

[41 

what is called the "aperture field method" for analyzing reflector antennas. 
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A REFLECTOR ANTENNA 


Fig. I 
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This approach is no less rigorous than the boundary value problem approach,' 
at least in principle. However, the aperture field method yields a great deal 
of insight into the operation of a reflector antenna. 

The aperture field method may be thought of as a solution to the boundary 
value problem with an infinite sum of partial solutions. One begins by calcu- 
lating the fields that would be produced by the feed in free space. However, 
the feed is not in free space. Its fields propagate to the reflector where they 
set up currents and secondary fields. These secondary fields in turn affect 
the feed, so the solution for its fields and currents must be modified. This 
effect is called feed reaction. The perturbation caused by the feed reaction 
then affects the reflector currents and fields. This iterative procedure could 
be continued until sufficient accuracy were attained. 

At first glance this procedure might seem more complicated than the 
general boundary value approach. However one can usually ignore the feed 
reaction in the case of a prime focus antenna. In addition, the feed reaction 
is of only minor importance in a Cassegrain fed antenna. 

The neglect of the feed reaction can be seen for a prime focus instru- 
ment from simple power considerations. Suppose the feed radiates power 

P as shown in Figure 1. The Pointing vector decreases according to the 
r 

2 

1/r law until the wave strikes the reflector. As will be shown later the 
reflector will columate this power into a plane wave. The power density 


D. 


of the reflected wave on the axis is thus 


— P where D. is the direc- 
2 r f 


4ttF 

tivity of the feed and F is the feed to reflector distance. The fraction of 

A f 2 

the radiated power absorbed by the feed will be approximately f J ■ 



8 


where is the effective area of the feed. For ordinary feed dimensions 

_ -4 

and distances F this ratio is less than 10 

Since the feed reaction is ordinarily so small the analysis of a large 
reflector antenna is a simple two step procedure. The first step is to 
analyze the feed as an antenna in feed space. The second step is to calcu- 
late the secondary fields from the reflector treating the feed fields as 
sources. The analysis and design of feeds is an art in itself, and has 
been discussed extensively^^’ Feeds will be given a minimum 

of attention in this paper. 

There have been two approaches to calculating the secondary fields 

r g] 

of the reflector. - One method, discussed by Silver 1 , calculates the 
secondary fields from the currents indiced in the reflector. In the other 
method, Silver , assumes the waves which are incident from the feed 
bounce off the reflector according to geometric optics. The far fields may 
then be calculated using an equivalence theorem. Both of these methods 
give essentially the same results. The latter viewpoint is the one that will 
be followed here. 

C. Geometric Optics 

The geometric optics approach is a simple way of determining the 
result of a wave striking a reflector. This method takes advantage of the 
fact that reflectors are typically hundreds of wavelengths across. The wave 
launched by the feed may be considered a scalar wave of the form^ 

§(r) = E(r) e“^ Y(r) . 


( 1 ) 
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The function § must satisfy the wave equation 

V 2 $ + k Z 2 = 0. (2) 

If one substitutes Equation (1) into Equation (2) he obtains the equation 

jVTj 2 =k 2 + -l JL • (3) 

In order to gain some insight into the implications of Equation (3), 
suppose the y and z derivations are zero. Then Equation (3) may be re- 
written 


2 


(4) 


From Equation (4) one can see that for k to dominate the right hand side 
of Equation (3) the per unit change in amplitude grading must be small on 
a wavelength scale. 

This condition is readily satisfied near the reflector of a large re- 
flector antenna. This obtains because the amplitude is determined by the 
feed gain function which is a broad slowly changing function. It typically 
changes by a factor of three from the center to the edge of the reflector, 

but this is over a distance of several hundred wavelengths. 

2 

From this discussion one can see that the k term in Equation (3) 
dominates in the region near the reflector. In this case the phase change 
of ¥ comes solely from the distance traveled. Thus the waves behave as 


plane waves in the vicinity of the reflector. Consequently the waves striking 
the reflector are reflected with their angle of incidence and reflection equal. 
One calls this kind of wave propagation geometric or ray optics. 
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D. The Fourier Transform Relations 

Geometric optics applies to the local regions near the reflector, but 
as the waves travel away from the antenna they must suffer diffraction. 

Thus, geometric optics applies only in the region near the reflector. 

Silver^" ^ makes use of an equivalence theorem to account for the 
diffraction phenomena that must occur for any electromagnetic antenna. 

The theorem requires one to draw a closed surface around the antenna. 

On this surface one places electric and magnetic charges and currents to 
account for the sources inside. In this case the sources are the fields 
produced by the antenna. The uniqueness theorem, then, tells us that the 
fields calculated from these "equivalent" sources will be the same as those 
calculated by analyzing the antenna as a boundary value problem. 

The surface for the equivalence theorem is traditionally drawn in a 
particular way. Part of it is a plane perpendicular to the axis of the antenna. 
This plane is called the "aperture plane. ” The surface is then closed around 


the back of the antenna. 


The equivalent sources that are placed on the surface are calculated 

from the fields that can be calculated by geometric optics. A slight problem 

arises in treating the fields near the edge of the reflector. In this region 

the amplitude of the fields are changing rapidly and hence do not obey geo- 

[ 12 ] 

metric optics. Silver's 1- solution to this problem is to let the geometric 
optics fields continue on to the edge of the antenna where they end on a line 
charge in the aperture plane. This convention simplifies the mathematical 
expressions because it makes the aperture plane electrically charge neutral. 
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In any case the error caused by edge diffraction is small. 

Under Silver's convention there are no charges or currents on the 

equivalent source surface except for the aperture plane. Thus the part 

of the surface outside the aperture plane can be ignored. In practice the 

electric and magnetic equivalent charges and currents are not calculated. 

This intermediate step is eliminated and the far zone electric field is 

written directly in terms of the electric field in the aperture. 

[141 

This relatively simple relationship may be written as 


E r .(a ) = -j 
ff — r J 


k 

4tt 


-jkr 


± r x [<i: z + £ r ) x J E(x) 


jkL(x) Jk -r 
e — e 


x 



(5) 


In this equation E (a ) is the far zone electric field at a field point in the 
it — r 

direction of the unit vector, a , k is the wavenumber which is Zrr /\, r is 

— r 

the distance to the far field point, a is the unit normal to the aperture 

z 

plane, E(x) is the magnitude of the electric field in the aperture plane, 
and kL(x) is the phase of the electric field in the aperture. These quantities 
are shown pictorially in Figure 2. Throughout this paper integrals will be 
assumed to be infinite over the whole plane, with E(x) falling to zero for x 
outside the reflector area. This convention has the advantage that the inte- 
gral in Equation (5) can be thought of as the Fourier inverse transform. In 
addition, linear polarization in the x^ direction will always be assumed. This 
is really no loss in generality since the analysis can be done twice if linear 
polarization cannot be assumed. Thus Equation (5) becomes 
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, -jkr 
„ , . _ . k e 

E ,(a ) = -j 

ff — r 4 tt r 


a x 
r 


(a + a ) x a j 
Z r x I J 




J E(x| 


( 6 ) 

For large antennas Equation (6) can be simplified since the vector ex- 
pression is very nearly -2a near the axis. Thus, the far electric field has 

—1 

the same polarization as the aperture fields near the axis. The vector dot 
product in Equation (6) can be expanded into a scalar dot product. Let us 
define the angular components in the x^ and x^ direction to be u^ and u^ with 

sin u^ = cos cp sin 9 
sin u^ = sin cp sin 0. 

Now the dot product becomes 

a - x = sin u x, + sin u„x„. 
r — 11 2 2 

The angles 0 and cp are shown in Figure 2. The angular components u^ and 
u^ are almost always so small that sin (u) can be'replaced by its argument. 
With these simplificatiqns the gain can be written as 

2 

. Ilk - <■ . jkL(x) jku . x , .2 

G( *) = | .J E M - - dx| , (7) 

o in 

x 

where P.^ is the antenna input power and is the characteristic impedance 
of free space. This equation states that the far zone gain is the square of 
the Fourier transform of the aperture electric field.’ Equation (7) can also 
be rewritten in the correlation integral form:. 


G <H> = 4^7?“ 

o m 


j A(i) • i d ,. 


( 8 ) 
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where 


A(t) = J E(x) E(x + t) e J k t L feE) - L (2£ + d£ . (9) 

X 

Equation (8) states that the antenna gain pattern is the Fourier trans- 
form of the field autocorrelation integral, given in (9), The vector t is the 
correlation distance between two points. Figure 2 shows the vectors x and 
_t in the aperture plane. 

Equations (8) and (9) can be rewritten in terms of the gain of the feed. 
The feed gain is usually better known than the electric field on the aperture. 
This form is 

G(u) = ~ J Mr) e — (10) 

T 

where 

A(t) = J T (x) T*(x +jt) e j k t L (x) - L (x + jr)} 
x 

where 


T(x) = 


1 G f (6 ' , cp') 
4tt (f + r 2 /4f) 2 


jkt(x) 

6 l 


and 


The angles 8 ' and cp are the angles of a ray from the feed to the point x 
in the aperture. The phase function -L(x) is the phase error caused by the 
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feed system, while L.{x) is the phase error caused by the reflector surface. 

Equations (7) and (10) are of central importance to this study. It is 
through these equations that the far field gain patterns will be related to 
phasing errors. Phase errors are expressed in the function L(x) since 
this function represents the phase of the electric field in the aperture 
plane. A paraboloid reflector has the property that it transforms a 
spherical wave from its focus to a plane wave. Thus, if the reflector is 
a perfect paraboloid and the phase surface of the feed is spherical, the 
phase error function, L(x), will be zero. However, any physical surface 
must suffer imperfections on some scale. These surface imperfections 
are the chief limitation on the performance of a reflector antenna. Various 
models for phase errors caused by surface imperfections will be presented 
in subsequent chapters. The implications of these models for the important 
antenna parameters such as gain, main beam efficiency, sidelobe level. 


and beamwidth will also be discussed. 



Chapter III 


A Critical Review of the Literature 

There are many articles in the literature related to the limitations 

of reflector antennas. This literature survey will include (A) a discussion 

of general antenna tolerance concepts, (B) a review of two papers dealing 

with specific phase errors, (C) a discussion and critique of the classical 

statistical tolerance model, and (D) a discussion of phase errors seen in 

practice. Section(E) shows a useful decomposition of errors. 

A. General Reflector Tolerance Theory 

[ 15 ] 

Bracewell showed that the aperture efficiency of a perfect re- 
flector is unity when the power from the feed is uniformly distributed 
over the aperture plane. In addition, he showed that any other distribution 
of energy produces an efficiency less than unity. Thus, the best feed 
pattern is one which illuminates the reflector uniformly and falls abruptly 
to zero outside the reflector area. In practice, of course, it is impossible 
to achieve such a pattern. Some of the feed power will inevitably miss the 
edge of the reflector. This power is called spillover. Thus, the feed de- 
signer must strike a compromise between illuminating the reflector uni- 

r 161 

formly and having too much spillover power . Ludwig had defined 

spillover and tapering efficiencies that separate these two effects for 

[ 17 ] 

comparison. In some receiving applications the spillover is very 

important because it contributes to the total noise figure as well as de- 
creasing the efficiency of the antennas. 
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The efficiency of a uniform phase aperture with the spillover, tapering, 
and blockage losses considered will be called the design efficiency in this 
paper. It is so called because it is fixed by the design of the feed system. 
The efficiency of an actual antenna system will fall short of the design 
efficiency because the reflector will not produce a true uniform phase front. 
The amount that the true efficiency falls below the design efficiency is called 
tolerance loss and may be expressed as a percentage loss or in decibels. 

B . The Effect of Specific Phase Errors 

A few authors have calculated antenna radiation patterns for specific 
phase error forms. However, most of the discussed phase error forms 

r 

are too restricted to be used in practice. 

[ 18 ] 

Dragone and Hogg have calculated antenna patterns for antennas 
with phase errors of the form; 

M 

L(x) V a cos (2 tt mr) (11) 

— <iTT L-> m 

m=0 

where 

r= |x|. 

This is a circularly symmetric phase error. Of course, any circularly 
symmetric function L>(x) may be expanded in such a series since it is a 
Fourier decomposition. 

They worked two examples which had only a single Fourier com- 
ponent. In both of these examples the sidelobe level was enhanced in 
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a single angular direction. The enhancement was also further from the 

main beam for larger m. 

[ 19 ] 

Silver has investigated circularly symmetric quadratic phase 
errors where 


L te> = 117(5) • 


The constant 6 is the phase error at the edge of the aperture in radians. 

[ 20 ] 

This work was based on an asymptotic expansion derived by Spencer 

J21] 

which was later corrected by Phgot . Silver showed that this kind of 
error causes the antenna beam to be broadened and raises the level of 
the near sidelobes. In Chapter III the quadratic type phase error will be 
considered in a much more general context. In particular it will be shown 
that moving the antenna feed will introduce an error much like the error 
given in the equation above. Thus, if an error like this were actually 
present in an antenna it could be substantially cancelled by merely moving 
the feed. 

C. The Classical Statistical Model 


Many of the papers on antenna tolerance theory are based on a 
statistical model for the phasing errors in antennas. This' model was 
first proposed by John Ruze in his doctoral dissertation. Ruze's model 
requires that one postulate a large number or ensemble of antennas. Each 
of these antennas is assumed to be identical except for "random" errors 
made in manufacturing the reflector. These manufacturing errors lead 
to a different phase error function L(x) for each antenna of the ensemble. 
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'Thus, L(x) is a two dimensional random process. In addition, Ruze 
assumed L(x) was wide sense stationary and Gaussian with a correlation 
function given by 

r l <s> = p *) 2 e-hl 2 /c 2 . 

The parameter e is interpreted to be the rms deviation of the reflector 
surface from a true paraboloid while c expresses the distance over which 
the phase errors are correlated. 

Since the function L(x) is a random process in this model the gain 

pattern, G(u) is a random process. This conclusion follows from the fact 

that integrals of random processes, if they exist, are random variables. 

The fact that G(u) is a random process poses a problem since the 

results of the theory must be applied to a single antenna. The solution 
F221 [23] . [24] 

adopted by Ruze 1 , Bao L , Shifrin , and others is to calculate 
the ensemble average radiation pattern. This is accomplished by taking 
the statistical expectation of the gain function. In this case one is aver- 
aging G(u) over all the antennas in the ensemble for each value of u. The 
average radiation pattern may be written from Equation (10) as 

E G(u) = £ J j E + T(x) Tfe + r) e'*- ' * dr. 

jr ££ 

( 12 ) 


where E stands for statistical expection. 
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In Equation (12) the expectation has been brought inside the integral. 
This is permissible since the integral of T(x) is bounded. Taking advantage 
of the stationarity of L(x) one can rewrite Equation (10) as 


EG(u) =~ 
— TT 


f( T ) A( T ) ^ d T 


(13a) 


where 


f( T ) = E e' 


jk[L(u) - L( t )] 


(13b) 


and 


A(r) = j T(x) T(x + r) dx. 


(13c) 


The function A(t) is the autocorrelation of the aperture field and is inde- 
pendent of any phasing errors. 

Another consequence of the stationarity of E(x) is that L(x) may be 
assumed to have zero mean. If it does not have zero mean then another 
random process say L/(x) may be defined as 


L'(x) = E(x) - m. 


where 


m = E L(x). 

Consequently one can easily see 


jk[L(o)-- L(t)] _ jk[L'(o) - L'( t )] 
e - e 


The function f(r) in Equation (13b) is a characteristic function. 
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A characteristic function of a random vector is in general 
f(w) = E e j - 2 

If one lets the y-vector be [L{o), L(t)] then the expression given by 
Equation (13b) is f(k, -k). Since L is by assumption a Gaussian random 
process with zero mean its characteristic function is given as 


f(w) = e 


- 1 / 2w Bw 


where B is a matrix given by b . = Ey.y.. Thus 

1 J 


B = (2e) 


1 e-UI 2 /^'’ 


.-111 /' , 


This is actually the definition of a Gaussian random process given in many 
[25] 

texts . Thus, by simple manipulation one has 


jk[L(o) - L(t)] _ -k [ R(o) - R(t)] 

-£_j © — 0 i 


(14) 


Substituting the form of R(t) into Equation (14), Equation (13) becomes 


E G(u) = 




- T 

i — 


2 , 2 

/ c 


A(t) e 


-j.k(u 


— d T . 


(15) 


Equation (15) is of fundamental importance since it is a completely 
general expression for the average gain pattern under theRuze statistical 
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model. Most authors who have dealt with this statistical model have 
derived their results under some condition which restricts the para- 
meters e and c or the angular variable u. 

1. Small Correlation Regions 

Ruze in his original work on the subject restricted himself 
to the case where c « D, where D is the diameter of the antenna. Ruze's 

r ? bi 

result L may be stated as 


E G(u) = G (u) e 
— o — 


~(o) 


_ , . .2 2 2n 

n=l 


(tt c p) 


n X 


(16) 


where; 

p - sin 9 

G^ = phase error free gain 

4n e 

CT = ~ ‘ 

The first term in Equation (15) has the same angular dependence as the 
phase error free or design pattern G . The second term is broader in 
angular extent. Its beamwidth is determined by the correlation distance, 
c. 


Equation (16) is interpreted to mean that energy is taken out of the 

main beam of the antenna and scattered into the sidelobe regions. Figure 

2 

3 is a normalized radiation pattern graph for c = 0. 05D with (4tt e/X) as 

[ 26 ] 

a plotting parameter. Figure 3 is borrowed from Ruze. As one can 

see from this figure, the first term in Equation (1 6) dominates on the axis 
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for e« X. Thus for axial gain considerations Equation (14) is often 
shortened to the form 


E-2i£)=e 
G (o) 
o 


/■4tt 

- f crv 


(17) 


Equation (17) is often used to extrapolate the measured peak gain 

» T 27 ~\ r 28T 

of particular antennas, to different wavelengths L ^ ’ . For this pur- 

pose Equation (17) is usually written in decibel form with the expectation 
being ignored, i. e. 


N_ = 680 (e/X) 2 , (17a) 

2 

where is the tolerance loss in decibels. If N^_ is plotted against (1/X) , 

a straight line should be formed if the tolerance loss follows Equation (17a). 

[291 

This graph has been called a Ruze diagram . The usual procedure is 
to plot several measured values of the tolerance loss on the Ruze diagram 
and fit a line to the data. Then the tolerance loss is estimated by the line. 

The procedure outlined above suffers from three shortcomings. 

First, it is a misapplication of Equation (17a). Since Equation (17a) applies 
to an ensemble average of antennas, the axial gain of a particular antenna 
cannot be expected to follow the law that Equation (17a) predicts, because 
the experimental outcomes of a statistical sampling scatter about the en- 
semble average. Secondly, Equation (17a) was derived under the assumption 
of a small correlation distance. As discussed in Chapter IV and pointed 
out by Ruze^" , errors often have large correlation distances in a practical 


case. 
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The third shortcoming relates to the mathematical interpretation 
of the Ruze model. The experimentally measured "gain" of an antenna 
is always the peak gain. Thus the measured peak gain should be com- 
pared with the ensemble average of the peak gain. The average of the 
peak gain is greater than the peak of the average gain function given in 
Equation {16). This happens because the peak gain may appear off the 
axis for some of the sample functions of G(u). The ensemble average 
peak gain should be written 

E sup {j j T(x) T (x + _r e jk ^ L( - ) ' L( - + -I )] e ~ jk - ' -dxdr) . 
u 

— _T X 

This effect will come into play most strongly for large c. This occurs 
because long correlation distance errors are most likely to change the 
pointing of the beam. 

In order to illustrate these three points an example has been con- 

2 

structed. The antenna for this example is a square with N independent 

regions as shown in Figure 4. Each region in this example may be thought 

■ [31] 

of as an individual reflector panel as discussed by Bao. The phase 

in each region is given by a Gaussian random variable, L.., from a random 
number generator. The formula for the tolerance loss of the peak gain as 
a function of frequency is 


v f > = -p 

x,y 


sin x sm 


j ~ f L.. ,2xm .2yn .. 
c n 1 .1 .2 


m=ln=l 


The variance of L„ is chosen equal to c/4n. Thus, the antenna gain maxi- 
mum (4. 3 dB loss) should occur when f = 1 according to the Ruze formulation. 
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Figure 5 shows six experimental outcomes or sample functions on 
a Ruze diagram. The antenna for this example has nine independent re- 
gions. It is felt that nine is a realistic number of regions for a practical 
antenna. As one can see the tolerance loss does not follow -any definite 
rule for large losses. In particular the losses deviate significantly from 
the line predicted by Equation (17a) in the large loss region. This is 
caused partly by leaving out the scattering term of Equation (16). How- 
ever, it is clear that the individual sample antennas not only fail to behave 
as Equation (17a) predicts but fail to follow etny definite law. 

Figures 6 and 7 are corresponding figures for antennas with 36 and 
100 independent regions respectively. It is clear that Equation (17a) gives 
better results as the number of independent regions on the aperture is 
increased, i.e. c decreased. In this case the gain integral, Equation (7), 
is actually an average of a large number of independent random variables. 
Thus, the variance of T\^(f) becomes small by the Central Limit Theorem. 

In this case a particular sample function T^(f) should agree with the ensemble 
average, E T^(f). 

It should be noted that small correlation distance phase errors do 
not guarantee that the antenna will follow the law given by Equation (17a). 
The assumption of Gaussian phase errors plays an important role. For 
example, if uniformly distributed phase errors are assumed instead of 
Gaussian errors, the ensemble average of the tolerance loss will increase 
faster than Equation (17a) predicts until the scattering term takes effect. 
Figure 8 shows six sample functions from an antenna with 100 correlation 
regions and uniformly distributed errors. 
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Fig. 5 




Tolerance Loss ( dB ) 
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Fig. 6 
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Fig. 7 



To! erQnce Loss ( dB ) 
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Fig. 8 
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2. Larger Errors and Correlation Regions 
L 32] 

Zucker derived a form for the axial gain similar to 

Equation (16) where he considered larger correlation radiuses. This 
result may be stated as 

2 = 


-a 


2n 


E G(o) = e 


•CT G (o) + e '2tt c"\ V o ; , , 'i , 1n ^ 

° V-rJ L STT L 1 ‘ A nJ (i8) 

n=l 


2c 


where all of the symbols have the same meaning as before and 


The parameter c must be such that 



» 2 


In the limiting case when c -* 0 this form agrees with the Ruze expression 
given in Equation (16). This form for the axial gain has not been widely 
used. 

[ 33 J 

Bao has turned to numerical calculations to see the effects 
of larger c on the radiation pattern. Some of the patterns which he calculated 
are shown in Figure 9. For small correlation distances compared to 

the antenna diameter Bao's results compare favorably with Ruze. For 
large correlation distance phase errors the sidelobes become very indis- 
tinct and merged with the main beam. For c comparable with D the relative 
power became quite high at the angular position of the first sidelobes in the 
design pattern. In addition there begins to be considerable beam broadening. 
This does not happen for short correlation distance errors unless the rms 
phase errors are so large that the main defraction beam disappears. 



RELATIVE POWER, dB RELATIVE POWER 
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o , C = 0/ 20 { i } e = 0.025 X 

b C = D/ 10 { (ii) e =0.050\ 

(Hi) € =0.075 X 
c C = D/5, e = 0. 075 X 

uniform 

— tapered 

6 is the r.m s surface error. 

C is the correlation interval. 

D is the aperture diameter. 

AVERAGE RADIATION PATTERNS 
OF CIRCULAR-APERTURE 
ANTENNAS WITH f (r) = (hr2) 
AND D = 360 X . 

Fig. 9 
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Scheffler has given an asymptotic form of Equation (16) 


fox large rms errors. This is 



where all of the symbols have the same meaning as before. The corre- 
lation radius and rms phase error determines the beamwidth here rather 
than the diameter to wavelength ratio. This simply means that the level 

of the main defraction pattern has sunk below the level of the scattered 

[35] 

energy. For even larger values of a Equation (19) become s[ 




( 20 ) 


Now both the peak gain and beamwidth are independent of wavelength and 
antenna diameter. The ratio e/c has been interpreted as the average slope 
of the phase error function. 

3. General Discussion of the Statistical Model 


Many of these results can be seen in an intuitive way by con- 
sidering Equation (15) as a Fourier transform. Specifically the autocorre- 
lation function A(t) and the design gain pattern are a Fourier transform pair. 
These are shown side by side at the top of Figure 10. 

The effect of the Ruze model is to multiply the autocorrelation 

I 

„ function by the function f(r) where 
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. .“|t| Z /c 2 


“O [ 1 


] 


f (r) = e 

This function is complicated in form, but it is almost equal 
to the simpler form 


f< T ) = a + (1 - a) e"l-l /c 


( 21 ) 


where 


2 

a = 

This simpler form has the same first three taylor series terms as well 
as the same asymptotic value as j t [ -» °°. Actually the form given implicitly 
for the decorrelation of L(x) by Equation (36) is just as valid as the form Ruze 
assumed. The advantage of this new form is that it may be easily Fourier 
transformed. Its transform is an impulse of heighth a plus a Gaussian of 
width /47tt7c and height (1 - a) . 

Four examples of the function I(_t) with various values of c, and e 
are given under A{jr). Their Fourier transforms are shown to their right 
under G^fu). The product of A(jr) f(jr) is the function actually transformed 
in Equation (13). However, multiplication in the correlation domain corres- 
ponds to convolution in the angular power spectrum domain. For this 
discussion it is easy to see the results of the convolution. 

The first example is the case of small a and small e. The trans- 


form is a relatively high impulse and a broad Gaussian under the impulse. 
The convolution of the error free pattern with the impulse merely replicates 
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and reduces the error free pattern. The convolution of the pattern with 
the broad Gaussian essentially replicates the broad Gaussian since the 
pattern is nearly an impulse with respect to it. The net result is a pattern 
like those shown in Figure 3. It has the main diffraction beam essentially 
unchanged and a broad skirt in the sidelobe region. This is the case con- 
sidered by Ruze.^ 

In the second example the a is large but c is still small. In 

this case the main diffraction beam sinks down below the scattering level 

[37] 

and pattern is determined entirely by c and e 

In the third example one has large a and c comparable with D. 
This is the case considered by Bao. The resulting pattern is like those of 
Figure 8. Here the Gaussian convolver is comparable in width to the main 
diffraction beam. Thus, it smears out the patterns and broadens them. 

The last example is that of large cr and c larger than D. In this 
example both the impulse and Gaussian beam are narrow with respect to 
the pattern so the error free pattern is merely replicated. This case is 
of only academic interest since a correlation distance much larger than the 
antenna means that there is no phase error. 

D. Phase Errors in Practice 

The relationship between the characteristic function, f(T), in the 
previous section and the resulting average antenna pattern is fundamental. 
This relationship may be stated from Equation (13) as 


G(u/\) =— F[A (t) f( T )-] 

— TT O— — 


( 22 ) 
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where F represents the Fourier transform. The fact that f(jr) was derived 
under a statistical model and G(u) was an average pattern is incidental. The 
function f(T) may be defined in a completely general way as 

J'j T H T(x + T) e-> kCL(x) - L(x + t)] dx 

X 

f(T) - — 

J T (x) T (x + r) dx 
x 

The pattern G(u) is now the observed pattern of the antenna rather than an 
average pattern. If there were no errors on the reflector then f(*r) would 
be unity. If any decorrelation takes place between points jr apart in the 
aperture then f(T) falls- below unity in magnitude. 

Equation (22) may be rewritten as 

k 2 

G(u/A) = — F [A ( T )-] * F[f ( t)] 

— TT O 

where * stands for convolution. Since f(T) represents the decorrelation 
of the fields on the aperture, one can classify the various kinds of decorre- 
lation from this equation and the pattern. 

The relationships shown in the previous section are still valid for 
this discussion if the pattern G(u) is interpreted as the observed antenna 
pattern instead of the average pattern. If the decorrelation of the fields 
takes place over a distance small with respect, to the diameter of the 
antenna, then the far zone pattern will have broad skirts. If the decorre- 
lation takes place over a distance comparable to the diameter of the 
antenna, the main beam will be broadened and the sidelobe structure 
generally less distinct. 
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This latter kind of antenna pattern is the one most often reported 

[ 38 ] 

in the evaluation of large reflector antennas. J. Schrarru reports 
that the beamwidths of the NRAO-36-foot antenna are about 10% different 
in the two principle planes. > This effect probably comes from a long 
correlation radius phase error smearing out the pattern in one of the 
planes . 

[39] [40] 

Cogdell and Bathker report almost identical phenomena 

in The University of Texas 16-foot antenna and the JPL 85-foot antenna. 
Both note that in one plane the pattern can be corrected by axially moving 

the focus. However, the pattern in the other plane is always broad. 

[41] 

Jacobs notes that at least one of the patterns of the Aerospace 

Corporation 15-foot antenna is always broad. He attributes the effect to 
the atmospheric scintillations. However, it should be noted that the 
broadness of the pattern could well arise from phase errors in the 
reflector. 

[42] 

Bathker gives a mechanically measured error contour map for 

the JPL 200 foot antenna. This map shows that the reflector is divided 
into six pie sections with each section having an error of the opposite 
sign from the two adjacent sections. This is clearly a long correlation 
distance error. 


Thus, many of the phasing errors usually seen in practice have the 


characteristic in common that they cause decorrelation over distances 
comparable to the diameter of the antenna. This effect seems to occur 



because the dominant error producing effects, such as thermal and 
gravity loading affect the antenna structure as a whole. 

E. A Natural Decomposition of Errors 

As discussed above long correlation distance errors scatter 
energy into the main beam and broaden it. This broadening can be 
used to separate the effects of long correlation distance errors from 
the short correlation distance "random" errors. 

It might at first seem artificial to classify errors as either of 
long correlation distance, comparable to antenna size, or of short corre- 
lation distance, much smaller than antenna size. One tends to exclude 
the intermediate case through a consideration of the sources of antenna 
errors. In the literature the following sources of error are identified: 

(1) gravity loading, (2) thermal distortion, (3) initial rigging or adjust- 
ment bias, caused for example by master template error or error in 
optical alignment instrument, and (4) manufacturing error or residuals. 

The first three sources affect the entire antenna structure and hence will 
produce errors which are correlated over the entire antenna. 

The fourth source will produce errors that tend to randomize over 
a small distance. In the case of a rigged antenna, in which the antenna 
figure is fixed by many mechanical adjustments, errors would decorrelate 
in a region controlled by a few adjustment bolts at most. If optical targets 
are placed on the surface to monitor adjustments, a similar situation would 
obtain. For a machined surface, errors would be determined by digitizing 



41 


residuals and deflections of the surface under the force of the machining 
tool. In all these instances, the correlation distance will be comparable 
to the micro structure of the antenna. 

The tolerance loss of an antenna is defined as the ratio of the 
true efficiency to the design efficiency, T)/T] . Let us break the tolerance 
loss into two factors, i. e. 


JL., ^ 

11 l s H 
o o 


11 


l . 


where T) is the tolerance loss from "random" errors and-—- is the tolerance 
3 \ 
loss from long correlation distance errors. The parameter TL is the effi- 

ciency of the antenna if only the long correlation distance phase errors 

were present. 

[ 43 ] 

If the main lobe is modeled by a Gaussian function then one can 
easily show that the efficiency is inversely proportional to the product of 
the beam widths. This proportionality holds because the integral of the 
pattern must be a constant. The appropriate efficiency for this discussion 
is T\ since it is assumed that the long correlation distance phase errors 
cause broadening of the beam, i. e. 


’‘"•A' 


If there were no phase errors, then T| should become 1] , so 


T] or 


o 0 *0 * ’ 


where 9^* and 0^* are the design beamwidths. 
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The expression for the tolerance loss can now be written 


n W 

U \ e e 

o 12 


The parameter T) is the loss from random errors, 
s 

[441 

given by the factor derived by Ruze , i. e. 


so it should be 


2 

^"4n e N , 

where e is now interpreted as the short correlation distance manufacturing 
error. Thus, one has 




_! 

T1 


= e 


-(&£) e *e » 
VW 12 


6 1 9 2 


Rewriting this expression in decibel form, one obtains 


H y 2 

N + 10 log ( i + log ( —=^) = 680 

t s ioV.9 ®iove^*y ’ 


(23) 


where N is the total tolerance loss in decibels, 
t 

2 

If one plots the left hand side of Equation (23) against f , then one 
should sense the losses caused by the short correlation distance manu- 
facturing errors. The graph which has been described is plotted in 
Figure 11 for The University of Texas 16-foot antenna. The value of e 
obtained is . 07 mm which compares favorably with the value of . 06 mm 


that the manufacturer quoted. This later value was measured relative 

to a template under ideal conditions when the antenna was being constructed. 




Chapter IV 

An Astigmatic Phase Error Model 

A. introduction 

The dominant phase errors in reflector antennas are often errors 

which have a large correlation distance. Direct surface measurement by 

optical or mechanical methods usually reveal only a few correction re- 

[ 45] , [461 

gions on the reflector ’ . In addition, antennas with significant 

[47 ] , [48 ] 

tolerance loss usually have broad beam widths and high sidelobe levels. 

In Chapter III it was pointed out that broad beamwidths and high sidelobes 

are evidence for large correlation radius errors. 

Often the phase errors in a reflector antenna system vary with the 

[ 49 ] , [ 50 ] 

antenna pointing . This effect is thought to arise from gravity. 

The gravity vector loads the antenna backup structure in different ways 

depending upon the orientation of the antenna structure. One might expect 

gravity loading to effect the structure as a whole and thus produce slow 

systematic errors. Such errors have been borne out by a structural 

[ 51 ] 

analysis done with a digital computer on the Haystack antenna. The 

calculated error was in the form of a four leaf clover with opposite leaves 
having the same sign error. The purpose of this chapter is to give a 
realistic model to this four-leaf-clover-type phase error, which is called 
astigmatism in optics. Several properties of the antenna pattern of an 
astigmatic reflector will be shown. In addition special emphasis will be 
given to determining the model parameters from experimental data. 

Both pattern range measurements and astronomical measurements will be 
considered. 
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Astigmatism is a well known phenomenon in optics. However, 

* 

different parameters are of interest in. antenna theory than in optics. 

Most antenna systems have only a single detector at the focus, as opposed 
to a photographic plate which is actually an array of detectors. An object 
must be observed in such a system by scanning the beam through it. Thus, 
the parameters of interest for an antenna system are the gain and the beam- 
width. In the case of an optical telescope one is interested in the extension 
of a star image caused by the astigmatism. 

B . The Model 

The model that will be considered is a model for phase errors in 
the aperture plane. This model may be stated as 


L r <*) 


2 2 



a 

mn 


n=0 m=0 


n m 



R m + n 
o 


(24) 


The symbol means a phase error caused by the reflector. The co- 
efficients a are model parameters and have units of length of phase 
error at the edge of the aperture. This expression may be thought of as 
a Taylor's series expansion of the actual phase error function. All the 
terms of less than third order are included. 


Only the second order coefficients of Equation (24) are significant 
to the far zone gain pattern. The coefficient a ^ does not appear in the 
expression for the far field radiation pattern while a^ and a^ represent 
linear phase gradings. As shown in Chapter II, the far zone electric field 
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is the Fourier transfor-m of the aperture electric field. Thus, the linear 

phase grading on the aperture represents translation of the beam in angle 

[52] . r . 

space. However, antenna properties such as peak gain, efficiency,. 

main beam efficiency, beamwidths, and sidelobe levels remain the same 

when the beam is translated. Thus translation of the beam is unimportant 

and will be ignored. 

C. Phase Errors Introduced by Moving the Feed 

The position of the feed is one of the fundamental parameters of a 
reflector antenna system. The feed position affects the phase front that 
will appear in the aperture plane. Thus, all of the important antenna 
parameters depend significantly on the feed position. 

The absolute position of the feed is seldom accurately known. More 
often, however, one has available a method of moving and measuring the 
feed position on an incremental basis. In the following an expression will 
be derived for the phase error caused by such an incremental displacement 
of the feed. 

Let the feed be moved in a coordinate system (A, B, C) as shown in 
Figure 12. The origin of this coordinate system is taken to be at the focus, 
but the final result would be unchanged so long as the origin were anywhere 
in the region of the focus. The length of a ray from the feed to the aperture 
plane is approximately 



1 

4f 


2 

r + 



A) + (x - B) + (F + C - 


4f 


2 2 
r ) . 


where 




Fig. 12 
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2 2 2 
r = X 1 + x 2 • ' 

For the discussion of this model lateral focus movement will be 
ignored, so A and B will be set to zero. However, it must be pointed 
out that lateral focus corrections should be an important part of any 
focusing procedure. Lateral focusing will be discussed as it relates 
to this model in. the next section. 

Assuming C « F and using the binomial expansion of the square 
root one can show that 


h ® = - c ! (S) +c 


_ “ 

1 C—') + 

4 V2F J 


where all of the constant terms have been neglected. The quadratic term 
is by far the dominant term here. For F/D = 0. 5 the quartic term is only 
1/4 of the quadratic term at the edge of the aperture. In the interior of 
the aperture where the field intensity is greater the quartic term vanishes 
comparatively. Thus, only the quadratic term will be retained in this dis- 
cussion. This assumption is consistent with the dropping of the third and 
higher order terms from the reflector error function. 

D. Total Phase Error 

The total phase error is. the sum of the error caused by moving the 
focus and the error arising from the reflector. This may be written 


L(x) = L f (x) + L^(x). 


Substituting the appropriate expressions for L^ and L^, one has 
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2i 2 

L( - ) = [ a 20 X 1 +a ll *1*2 +a 02 X E +C 


5 i 2 ^ 2 ^/t> 2 

T <IC 1 + x 2 »j /E o ■ 


64{F/D) 


The constants a , a , and G are redundant so it will prove convenient 
20 02 ± 

to eliminate one of them. Let the new set of constants be ct, p, T where 


L(x) = or +p x i x 2 ~ 01 x 2 + r K ( x j + x 2 ^ R 0 ( 25 ) 


where 


tr = ■ 


a 20 " a 02 


r = c + 


a 20 + a 02 
2 K 


P =a ir 

5 2 

and K= 64VF^ • 

The constants o’ and p depend only upon reflector errors while T is directly 
proportional to the focus position. T also depends on the reflector errors. 
However since the reflector errors do not change T can be thought of as 
the focus position variable. 

[53] 

An interesting special case of Equation (25), given by Silver 1 * , has 
been discussed in - Chapter III. He considered a quadratic reflector error of 
the form 


L(x) - Cr 2 

where C is a constant. 

Thus in Equation (24) a^-, = a 2 Q an< ^ a ]q = 0- The phase error given by 
Equation (25) can now be made identically zero by moving the feed so that 
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L = 0. Thus, since this kind of phase error can be corrected by simply 
moving the focus position, the quadratic reflector error given by Silver 
is degenerate. Of course, it should be kept in mind that moving the feed 
introduces higher than second order phase terms. 

E. Model Properties 

As discussed in Chapter III phase errors most often cause decorre- 
lation over distances comparable to the diameter of the antenna. Hence, 
the astigmatic phase error model which has been discussed has a certain 
intuitive appeal. In this section several model properties concerning the 
far zone radiation pattern will be given. Through these properties, one 
can determine whether this model is a reasonable representation of the 
phase error function. 

In order to state some of these properties in their most general form, 
the concept of an antenna pattern contour map must be introduced. A contour 
map is a two dimensional representation of the radiation pattern of an 
antenna. It has two rectangular coordinates u^ and u^ which correspond 
to angular directions u^ and u^. Lines are drawn in these coordinates 
which describe the locus of a specific pattern level. 

For this discussion it will prove convenient to define two symmetry 
concepts. One is symmetry about a line through the origin at an angle cp. 

The angle cp is an angle measured counterclockwise from the u^ coordinate. 
Symmetry about this line means that the mirror image of a locus on one 
side of the line appears on the other side of the line. Another concept that 
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must be defined is symmetry of the map about the origin. This will mean 
that all of the contour loci are symmetric about the origin and that the gain 
function is even, i. e. 

G(u) = G(-u). 

These five properties may be thought of as symmetry properties of 
the far zone pattern. Some symmetries can be seen at any focus position 
while others require focus movement. The proofs of these properties are 
given in Appendix A. 

The first two properties are symmetries which can be seen at any feed 
position. The first of these may be stated as follows: 

Property 1. If the illumination function E(x) is symmetric about the origin, 
i. e. E(x) = E(-x), then any contour map is symmetric about the origin for all 
feed positions and reflector errors. That is 

G(_u) = G(-u), for all a, [3, and F. 

■ This property also implies that a pattern cut taken in any direction will be 
symmetric. 

The condition E{x) = E(-x) is a very mild one. It is difficult to imagine 
any feed system for which the illumination function would not be even, since 
most feed systems are symmetrically constructed. 

Figure 13 is an example of a contour map of the pattern of a reflector 
with astigmatism. This figure illustrates property 1 holds in the finest 
detail. Figure 13 and several figures that follow in this chapter were cal- 
culated in a digital computer with a program that will be discussed in 



52 - 


Appendix A. The illumination function is the illumination of a horn feed 
of dimensions 1,25 1 x 0.9 X on a paraboloid with an F/D = 0.5. 

The converse of this simple property is perhaps more useful, since 
if the contour map is not symmetric the model does not hold. In an actual 
antenna evaluation, of course, one does not expect property 1 to hold with 
arbitrary accuracy. 

Pattern asymmetries can sometimes be corrected by moving the feed 
laterally. It is easy to show that the pattern assymmetry is caused by an 
odd phase error function, that is a phase error such that 

L(x) = -L(-x). 

Moving the feed laterally produces an .incremental phase error that also 

has the odd property. This incremental phase error can sometimes be 

used to cancel the odd error in the reflector. The. focus position can be 

moved laterally until the level of the first two sidelobes are equal in each 

[54] 

of the two principle planes . It is easy to show that the feed should be 

moved in the direction of the lower sidelobe. For example, if the lower 
sidelobe is East of the main beam, then the feed should be moved East. 

If the patterns are still asymmetric after the sidelobe levels .are 
made equal then there is a large odd phase error. In this case the 
astigmatic phase error model is clearly inappropriate. 

The second property shows that the contour map should be symmetric 
about two perpendicular lines. Property 2 may be stated as follows: 
Property 2. If the illumination function E(x) is only a function of r, the 
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the radius from the center of the aperture, then a contour map of the 
pattern is symmetric about the lines 

tan 2cp = ~ 

2a 

for all feed positions. 

Property 2 gives a further fact that can be used to determine 
whether a given reflector has astigmatism. If the property is satisfied 
then the relationship between a and 8 is easy to establish. One need only 
plot a contour map of the antenna and draw the cp direction on it. Then a 
linear relationship between a and 8 is established by the equation given 
above. 

Figures 13 and 14 illustrate this property. In the first case the 
cp directions are aligned with the axes. In Figure 14, however, the cp 
direction is inclined to the axis by about -9°. As one can see property 2 
is not rigorously satisfied in Figure 14. This occurs because the illumi- 
nation function is not quite circularly symmetric. The illumination, 
however, 'is chosen to be realistic. It is the illumination produced by 
the horn feed discussed above. 

Property 3 shows that the radiation pattern produced by a phase 
error of o’ and 8 for some F is the same radiation pattern that would be 
produced by -o', -(3, and -T. Property 3 may be stated as follows: 

Property 3. If the illumination function E(x) is symmetric about the origin, 
i. e. E(x) = E(~x), then the gain is an even function of <y, B, and F taken 
together. That is 
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G(u r u 2> O', p, r) = G(u 1 , u , -O’, -p, -r). 

Figures 13 and 15 illustrate this property. These two contour maps are 
identical while Figure 13 was calculated with the negative of the parameters 
used in Figure 15. 

In a subsequent section a technique will be presented that predicts 

k 

the model parameters by matching the predicted radiation pattern of this 
model to the observed pattern. Property 3 makes the results of any such 
technique ambiguous. Specifically, it is impossible by this method alone 
to determine whether a certain reflector error a and g or its negative -or 
and -g is present in the reflector. 

Property 4 states that changing the sign of the reflector errors o' and g 
exchanges the role of the patterns taken in the cp direction. The u^ pattern 
becomes the u^ pattern and the u^ pattern becomes the u^ pattern. Property 
4 may be stated as follows: 


Property 4. If the illumination function E(x) is only a function of r, then 
the angular coordinates which are aligned with the cp direction, tan 2cp ~ > 

change roles when a and (3 change sign. That is 


where 


G(ui', u 2 7 , a, g, r) = G(u 2 ', u x ', 


-a, -8, D 



This property will be useful in a subsequent section. 
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Property 5 gives the most detectable consequence of primary 
reflector astigmatism. It may be stated as follows: 

Property 5. If the illumination function E(x) is only a function of r, then 

the angular coordinates which are aligned with the p direction, tan 2cp = , 

cot 

change roles when r changes sign. That is 


G(Ui', u^', of, p, T) = G(u 2 ', 


, a, -p, -r) 


where 


/ 

U 1 


cos cp 

- sin cp 

1 d 
N 

1 . ' 


sin cp 

cos cp 

— 


One immediately correlary of property 5 is that the axis gain is an 
even function of T l\. The axis gain is the case where u^ = u^ = 0. Figure 
16 shows the axis gain plotted against T with o' as a plotting parameter (p =0). 
This graph shows that the axis gain suffers little when a is small, but deter- 
iorates rapidly as a is increased. The axial gain is also much less sensitive 
to the feed position when of is large. 

Property 5 also requires the radiation patterns in the cp directions to 
exchange roles when F changes sign. This is perhaps the most striking quality 
of the model. As the feed is moved from -F to T, the u^ ' coordinate should 
become the u^ / coordinate, and the u^ ' pattern should become the u^ 7 
coordinate. Figures 15 and 17 illustrate this behavior. The illumination 
function for these maps was calculated for a horn feed as mentioned above, 
so the illumination is not quite circularly symmetric. That property 5 
is not satisfied exactly can be seen in the 30 dB contour. 




r '\ 


Fig. 16 
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Suppose the u ' pattern is broad while the u ' pattern is narrow 

-L Li 

at some focus position. Then there should be another focus position 

where the u ' pattern is broad while the u ' pattern is narrow. Figures 

18 and 19 illustrate this fact. Figure 18 shows the principle cuts of the 

pattern with ko = 2, 8=0, and V/\ = +0. 2, while Figure 19 shows the 

pattern cuts with ky = 2, 8 = 0, and T/X = -0. 2. The u^ pattern of Figure 

18 is the u pattern of Figure 19 and vice versa. 

This kind of pattern behavior has been noted by several authors. 

[55] 

Bathker in his evaluation of the JPL 85-foot Ground Antenna states, 


The hyperboloid was set 0. 500 in. toward the apex which approxi- 
mates the high elevation angle focus for maximum gain. Figures 
29 and 30 show the resultant patterns. The azimuth cut is seen to 
have a narrower main beam while the elevation cut is broader. . . . 
Although a corresponding hyperboloid position toward the vertex 
f was not tried it is suspected the inverse would be true, i. e. the 
elevation would be nearer an optimum focus and the azimuth plane 
would be defocused. . . . 

Bathker later confirmed that the paraboloid had a surface error that would 
fit the astigmatic phase error of Equation (25) very well by making mechanical 


measurements of the reflector figure. 

Bathker is not alone in reporting this effect. 


Jacobs and King 


[56] 


in reporting on the characteristics of the Aerospace 15-foot antenna state: 

. . . On rare occasions only moderate beamwidth broadening in the two 
orthogonal patterns occurred. However, on most occasions there was 
modest broadening in one plane and considerable broadening in the 
orthogonal plane. Sometimes this large amount of broadening occurred 
in the declination cut, and on other occasions it occurred in the hour- 
angle cut. On no occasion did extensive broadening occur in both 
planes during one measurement period. . . . 
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Fig. 18 
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Fig. 19 
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Jacobs and King interpreted this effect as arising from the atmosphere. 
However, it seems strange that the atmosphere could somehow scatter 
energy in the declination direction and not in the polar direction. Instead, 
this effect could quite possibly have come from slight movement of the 
secondary reflector and astigmatism in the primary reflector. 

F. Measurement of Model Parameters 

If adjustments are to be made based on this phase error model it 
is important to know the model parameters. Three methods will be given 
for determining these parameters. The first method depends upon the use 
of a pattern range facility to make a contour map and beamwidth measure- 
ments. The second is a method for measuring beamwidths astronomically 
by observing solar limb crossings. The third technique is a digital com- 
puter algorithm. It matches a predicted map of the antenna pattern to an 
observed map. 

1 . Pattern Range Method 

The model parameters can easily be measured if a pattern range 
is available. The first step is finding the line of symmetry on the contour 
map. The line of symmetry can usually be found from a single pattern 
contour. A 10 dB level has been found to be convenient, because it is sensitive 
to the cp direction yet it is not usually complicated by sidelobes. If one 
identifies the symmetry directions on the contour map, then a linear 
relation between a and f3 is established by Property 2. This is 


- 2ct tan 2cp = p . 


( 26 ) 
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It is possible that the contour map will be circular. However, at 
another axial focus position the eccentricity of the contour map should 
appear. If the contour map were circular at each axial focus position 
no astigmatism is present in the antenna. 

The parameter cv can be determined by observing the beamwidth 
change in the cp direction with the axial feed position. Figure 20 is a graph 
of the beamwidth in the cp direction vs focus position. The abscissa scale 
is T/\. The ordinate scale is the fractional widening of the beam over 
the design pattern beamwidth. The plotting parameter is key*, where k 
is the wavenumber. The relation between a' , o' and cp is given by 

of = a ' cos 2cp. (27) 

The parameter oi' is the value o' takes on when the (x^,x^) - coordinates 
are rotated so that the center term in the model. Equation (25), is elimi- 
nated. 

The value of a ' can be determined by plotting the observed beam- 
widths in the cp direction vs focus position on Figure 20. The parameter 
a' is then estimated by eye. This task is relatively simple because for 
phase errors greater than one radian at the edge of the aperture the beam- 
widths are very sensitive to axial focus position changes. 

If the measured beamwidths match the predicted beamwidths but with 
the coordinate labels reversed, then the sign of a 7 is negative. This fact 
follows from Property 3. 

With a' and cp determined, the parameters a and p can be calculated 
from Equation (27) and (26). 



BEAMWIDTHS 


Predicted Beam widens for 
F/D = 0.5 ond Horn feed . 


1 .30 


1.20 


1.10 
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The parameters a and 0 were estimated by this method for 
The University of Texas 16-foot antenna as discussed in Chapter VI. A 
contour map was plotted and the cp direction was taken to be - 18° . The 
cp direction beamwidths were then plotted, and the value of a' was esti- 
mated as 1.5. The values Equations (27) and (26) give for a and 6 are 
a = 1.21 and 0 = 1.75. This compares favorably with a = 1.13 and 0 = 0.73, 
given by the more objective computer minimization method. The discrepancy 
in 0 can be explained by the fact that the angle cp was clearly taken too large, 
see Figure 26, Chapter VI. 

2. Astronomical Method 

The method given above depends upon measuring the beamwidth 

change as a function of focus position in order to estimate the magnitude 

of the astigmatism. Beamwidths can also be measured astronomically if 

a pattern range is not available by observing solar limb crossings. This, 

[57] 

method has been used by Jacobs and King with good results. Also a 
solar limb crossing experiment is reported in Chapter VI. The beamwidths 
that were determined from this experiment compare favorably with the 
beamwidths measured on the pattern range. 

The passing of the solar limb through the antenna beam produces 
a record proportional to the integral of the antenna pattern. Thus, is is 
expedient to choose a model for the antenna beam. Almost any model 
would suffice since one is interested in only the gross pattern features, 
so a Gaussian pattern will be chosen for analytical convenience. In the 
(x^x^-coordinates of Figure 21, the antenna pattern model may be written, 
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G <*i' x 2> = °o e 


x- - Kt x„ 

-{(•VKrj)- 


(28) 


where G is the peak gain, 0 and 0 are the semi- beam widths, K is the 
o 12 

solar rate, and t is time. 

The sun is a thermal source, so the antenna smoothing 
equation applies. The antenna response is thus, 


CO 

J , 'j G(x i: X 2 >T(x r X 2 )d5C l dX 2' <29) 

— CO 

where T is the antenna temperature, and T(x.,x ) is the equivalent 

3 . * C* 

temperature at the point (x , x ). Substituting the form given in Equation 

* 6 

(28) for G into Equation (29) and assuming the sunhas a linear edge, one 
obtains 


T 

a 


= ~(t + T ") + ~( t : 

2 V sky sun/ 2 V sun 


- T ) erf 
sky/ 


f (Kt - cp) 


0 1 'V-Ln2 


(30) 


T ^ is the equivalent temperature of the sky, is the equivalent tem- 

perature of the sun, and erf(x) is the error function. Equation (30) states 
that the antenna temperature increases from the sky temperature to the 
solar temperature as the antenna beam crossed the limb of the sun. The 
scale of this transition is proportional to the half power beamwidth in the 
direction that the. antenna is driven across the solar limb. 

The procedure for measuring 8 is to let the antenna tempera- 
ture during the limb crossing be the data set, i. e. * 
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T (i) , i = 1 , 2, ... n. 

EL 

Then one estimates 0^ by minimizing the rms error between the data 
and Equation (30) with a digital computer. 

The experiment that has been outlined gives the equivalent 
beamwidth in the direction the antenna crosses the solar limb. Thus, 
in order to obtain a contour map one must take the beamwidth in a 
number of directions. This will determine the tp direction just as before. 
One can then determine the parameters by plotting the beamwidths in 
the cp direction vs axial focus position. 

This method is much more tedious than the pattern range 
method and requires the use of a digital computer.' However, it does 
not require the availability of a pattern range. This is important because 
some antennas do not have a pattern range. 

3. Computer Algorithm 

The computer algorithm method of estimating the model 
parameters depends upon matching a predicted pattern map with the 
measured radiation pattern of the antenna. The radiation pattern of a 
reflector antenna is sensitive to long correlation radius phase errors in 
the angular region near the main beam. Thus, even with only a small 
dynamic range there should be little difficulty in obtaining a relatively 
noise free measurement of any' astigmatism present. 

This algorithm calculates the antenna radiation pattern using 
the Fast Fourier Transform method presented in Appendix A. The pre- 
dicted radiation pattern can be computed to slide rule accuracy in less 
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than one second on the CDG 6600 computer. This high speed operation 
makes the iterative prediction of the model parameters possible since 
the predicted radiation pattern must be calculated hundreds of times 
in such an iterative method. 

The data for this experiment is a map of the antenna pattern 
near the main beam digitized at an N x N grid of points. The map is 
normalized to unity power level on the axis. This data set may be 
expressed as 


D(i, j) = G ' (u. , u.) 

^ J 

where 


u = [i-n/2 + 1] Au^ 

u, = [j-n/2 + 1] . 

J c 


(31) 


G ' is the measured normalized pattern of the antenna in dB units. The 
normalized axis gain G'(0, 0) appears one point to the right and above the 
center of the array. The sampling intervals Au^ and'Au^ should be chosen 
smaller than the maximum dictated by the sampling theorem 1- , i. e. 


iu r lu 2 



radians. 


The sampling intervals in the two directions need not be equal. 

The algorithm minimizes the objective function, F, over the 
four parameters a, b, a , (3, T where 


n n 


f = Z I [D(i ’ j} - 

j=l i=l 



The predicted antenna pattern, P(i,j), is calculated by the program of 
Appendix A. P{i, j) is an approximation to the integral 
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P(i, j) | J T<x r x 2 ) e jkL( - } e jk - * - dx j , 


x 

where u. ( and u, are defined by Equation (31). 
The phase error function L(x) is given by 


2 2 2 2 2 
L(x) = axj + bx^ + {a x 1 + g x^^ " a x 2 + T K(xj + x„, )}/R q 

The linear phase grading terms are included so that the predicted pattern 

can be translated to agree with the observed pattern as much as possible. 

To test the method and the program, the data array was set 

to the pattern of a uniformly illuminated constant phase aperture. This 

[59] 

data pattern was calculated analytically from the form given in Silver 
Then all of the phase error parameters were each initialized at 2/k where 
k is the wavenumber. All five parameters converged toward zero as 
expected. When the program ran out of time each parameter was less 
than 0. 01 /k. 

The program was then tried with real data. The pattern map 
was taken by J. R. Cogdell on January 30, 1969. The results of this is 
shown in Figures 22 through 24. Figure 22 shows the predicted and ob- 
served radiation patterns on the same scale. Figure 23 shows the observed 
radiation contour map while Figure 24 shows the best fit map. The observed 
map shows a considerable amount of beam squint which cannot be taken 
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into account by this model. However, this kind of error can sometimes 
be corrected by simply moving the feed laterally. 

It should be remembered, of course, that the results of this 
minimization procedure are ambiguous. Property 3 guarantees that the 
predicted pattern for a, £, and T will be the same as the predicted pattern 
for -u, -£}, and -F. This ambiguity can be resolved by plotting the beam- 
widths against feed position as shown in the previous section. 

Both a conjugant gradient iterative method and a simplex 
method were tried with almost identical results. The conjugant gradient 
method, however, was somewhat faster. 

The computer algorithm method for determining the model 
parameters is the most objective of the> three methods and probably is 
the best. However it is also the most time consuming because of the 
programming effort required. In addition, it requires both a pattern 
range facility and a large digital computer. 

G. Conclusions 

An astigmatic phase error model has been presented. It has been 
shown that the axial feed position plays an important role in determining 
the character of the total phase error at the aperture plane. Five model 
properties have been presented with examples. Property 5 is particularly 
interesting since it has been seen in at least three large reflector antennas. 

If astigmatism is identified in the reflector then one is interested in measuring 
it quantitatively. Three methods are given for doing this. The first and 
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third methods depend upon pattern range measurements. The former is 
a field method; the later uses a computer minimization technique. The 
second method is an astronomical method which can be used if no pattern 
range is available. 



Chapter V 


Reflector Efficiency Evaluation By Frequency Scaling* 

A. Introduction 

The useful bandwidth of a parabolic reflector antenna is limited by 
the precision of its surface and the integrity of its supporting structure. 
Surface deviations from a true parabolic degrade performance as the 
frequency of operation is raised [60]. The number of antennas in which 
such effects have been observed over a wide range of frequencies is con- 
tinually increasing. 

Most applications involving reflector type antennas require only a 
knowledge of the antenna properties as a function of frequency. However, 
there is an increased class of problems where precise knowledge of 
antenna characteristics is required. Space communications and radio 
astronomy provide significant examples where large antennas are used at 
many frequencies. The Ruze tolerance theory [61] is currently used to 
predict antenna efficiency at frequencies where performance has not been 
evaluated by direct methods. 


*A large portion of this chapter is taken from a paper of the same title 
submitted to the IEEE Transactions on Antennas and Propagation. This 
paper was co-authored by J. R. Cogdell and the author. 
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The Ruze tolerance model has successfully described the variation 
of small tolerance losses with frequency in several antenna evaluation 
programs [62], [63], The wide success of this model is somewhat sur- 
prising in view of the strict statistical assumptions which underlie the 
theory. As duscussed in Chapter III, the model assumes that the phase 
fronts in the antenna aperture may be described by a two dimensional 
Gaussian random process. It also requires that the random process be 
stationary and de correlate over distances larger than c, where c is much 
less than the diameter of the antenna [64]. 

As discussed in Chapter III, a number of large reflector antennas 
have been subjected to extensive mechanical error evaluations [65], [66]. 

The results are usually displayed in the form of error contour maps. Exam- 
ination of these maps confirm the impression that errors tend to be corre- 
lated over large regions of the antenna, with the result that only a few 
independent correlation regions exist in the antenna surface. There are 

further reasons to suspect that errors should correlate over large fractions 

/ 

of the antenna surface. Error producing effects such as gravitational 
loading and thermal expansion act on the structure as a whole, producing 
systematic deformation in the structure. As shown in Chapter II, the 
tolerance loss of an antenna will not follow any precise law if there are 
only a few correlation regions on the aperture. 

The purpose of this chapter is to describe the tolerance loss in reflec- 
tor antennas as a function of frequency without the statistical assumptions. 
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The theory is based upon measurements and estimates of peak to peak 
error. The major points of the theory may be summarized as follows: 

1. An efficiency measurement and error estimate at one frequency 

implies firm bounds on the efficiency at other frequencies. 

These bounds can be used with confidence because nothing 
about the statistics of the phase errors is assumed. 

2. Measurement errors scale down with frequency. Thus, a 

moderately accurate measurement of the fficiency at a high 
frequency can be used to yield a highly accurate inferred 
efficiency at a lower frequency. This approach has the added 
advantage that gain measurements at higher frequencies are 
less subject to multi-path interference. 

In the following section tolerance loss is discussed generally. A 
series form for the telerance loss is presented and discussed with reference 
to the Ruze theory. In the fourth section the theory of efficiency measure- 
ment by frequency scaling is presented in the form of two theorems. The 
first theorem gives a method for predicting the efficiency at one frequency 
given a measurement at another frequency with an electrically scaled-feed 
system. The second theorem extends the first theorem to include gain 
measurements at two frequencies, which are used to predict the efficiency 
at a third frequency. These two theorems are rigorously true only for 
prime-focus instruments with electrically scaled feed systems. 
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This assumption is required to make the illumination of the reflector 
the same at all the frequencies. The first theorem, however, is extended 
by an approximate method to cases where small changes in the illumination 
takes place from frequency to frequency. Cassegrain antennas are also 
considered by this approximate method. The use of these methods will be 
illustrated by published data. 

B. Tolerance Loss 

As shown in Chapter II, the efficiency of a reflector antenna, at a 
frequency f, may be expressed as 

Tt(f) = F J B(t) dr 

_T 

where 

B(t) = J E(x) E(x+t) + 1)1 dx 

x 

and 

F - — 

2Z P A 
o r g 

The function E(x) is the. magnitude of the electric field at point x in the 
aperture plane and is assumed known. The phase length function, L(x), is 
usually unknown and contains the phase deviations caused by reflector 
defects. B(^) is the autocorrelation function of the fields in the aperture 
plane. The other symbols are r|(f), the aperture efficiency as a function 
of frequency; k = 2nf/c; Tj = ohmic efficiency of the antenna, Z^J[X Q t e o * ; 
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P the radiated power; and A , the geometric antenna area. For .the rigor- 

g 

pus results it is assumed that the antenna is illuminated the same at all fre- 
quencies, such that E(x) does not depend on frequency. 

It is well known that the autocorrelation function C(t) is Hermitian, 

i. e. , 

C( T ) = C* (- t) 

Thus, when the integral over t is performed the imaginary part of the 
integral will vanish, and we may work simply with the real part. Hence, 

rj(f) = F J J E{x) E(x + t) cos k[L(x) - L(x+_r)] dx dr_. 
jr x 

One can define the sequency of coefficients A as 


2n 


.n Zn 

= F JJ E( - } { c [ u - ] " } - '<% <*T . 


TX 

such that one obtains an infinite series in even powers of frequency: 

CO 


p{f) = £ f 2n 
n=0 



{ 32 ) 


In this expansion the first term is clearly the phase error free efficiency, 
which can be calculated without knowledge of reflector errors. The 
second term will describe a loss in efficiency which is quadratic in fre- 
quency and will show the beginning effects of tolerance loss. The quartic 
term will be important as tolerance loss gets substantial and so forth. 

i 
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Other aspects of Equation (32) should be noted. The first is that the 
series can be terminated after a few terms and the remainder bounded. In 
cases where peak phase errors are not unreasonably large the higher order 
terms are of rapidly decreasing importance. With one or two measure- 
ments one can estimate the important coefficients and bound the efficiency 
characteristics in the useful region of the antenna. This procedure will be 
quantified and illustrated in the final section of this paper. Spencer 
gave a form equivalent to the first two terms of this series. He made no 
effort to bound the error of approximation, and did not deal with the effects 
of measurement errors. Equation (32) can also be related to the statistical 

ri on 

results of Ruze L . The Ruze theory is applicable to the average of 
Equation (32) over a large number of similar antennas. The effect of the 
statistical assumptions is to determine all of the coefficients of the aver- 
aged equation in terms of two parameters, the rms surface deviation and 
the decorrelation distance. 

The Ruze tolerance loss factor is of the same form as Equation (32) 
when expanded in a power series in frequency. The first effects of toler- 
ance loss, which are always quadratic in frequency by Equation (32), can 


therefore always be fit by the Ruze factor, e v ^ . This might explain 

the fact that the Ruze factor is useful in the small tolerance loss region, 
say less than 6 dB tolerance loss, for some value of g, which is often 
called the "electrical tolerance. " It is relevant to consider that the 


"electrical tolerance" is often significantly different from the measured 
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mechanical tolerance ’ 5 , while the Ruze theory requires that 

they be the same. Thus, our conclusion is that the wide success of the 
Ruze tolerance loss equation is due more to the functional form of the equa- 
tion than to the validity of the model. 

C. Notation 

Before stating the theorems, we shall define the notation used in the 

two theorems. These quantities are as follows: 

T1 = design efficiency {the efficiency with no phase errors) 

T|(f ) = true overall efficiency at {(unknown) 

r]'{f) = measured overall efficiency at f 

B. = peak measurement error bar at f. 
i i 

K(f) = peak-to-peak phase error at frequency f 

f = measurement frequency 

m 

f = first measuring frequency 

f = second measurement frequency 

f = frequency of estimation (frequency for which 

efficiency estimate is desired) 

S(f, x) = phase error function = kL(x) 

D. Theorems 

Theorem 1 gives a method of predicting the gain at frequency f, based 
on the gain measured at another frequency, f , and an estimate of the peak- 
to-peak phase errors. The error made in estimating the gain is from two 
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sources. One is the error in the measurement at frequency £ ; the other 

is the mathematical error made in the estimation. Theorem 1 may he 
stated as follows: 

Theorem 1 

If the feed system of a prime focus reflector antenna at frequency 

f is a scale model of the feed system at f , and 
e m 

| rj* <f ) - r)(f ) I < B 
1 m m 1 — 

and 

I 5 (f , x) - 6jf , y) I <K(£ ), all x and y 

then the relations 


T)\ 


<y < \ - P[T1 0 - n' (f m n (-j 5 -) 2 + T| o sup S+ P b(-j 2 -) 2 , for all P 

m 


e _\2 
m' 


and 


1 <f„) > % - P[ r, o - (f-J + „ o mf S - P E Qf- y, for all P 

m - m 


hold true. The set S is defined as 


S - {z:z = U x) - P (jr-y U (x); 0 < x < K(f m )} 
m . m 


The function U(x) is given by 


U(x) = cos (x) - 1 

The proofs of the theorems will be given in Appendix C. 

In the bounds for the efficiency given above, the terms n Q sup S and 

p inf S bound the mathematical error caused by approximating cosine 
o 
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while the terms + 



result from the measurement error. 


4 

The quantity 


m 

P is a free parameter, so the method is most effective if the bounds are 
appropriately minimized or maximized over P, which is normally different 
for the upper and lower bounds. 

In the case where f < f very good results may be obtained using 

e rn 

P = 1 for the upper bound and 


D Of K < £ m >) 

m 

U(K(f )) 
m 



for the lower bound. In that case one has the simple expressions 


and 


° (r K(£ m>) 

U(K(f ))■' £ % ~ T l ,(f m ) ' ^ 
m 


’I'V a "o ‘ C'r ) 2 £% ' + 


Theorem 1 may be applied to predict the efficiency of a prime-focus 
reflector antenna with a scaled feed system at either a higher or lower 
frequency, but it is most useful in the latter case. This is true since the 


factor 



scales down the measurement error in the case where f 

e 


m 

is less than f . The University of Texas 16 -foot antenna is a good example. 
The required data is as follows: Cogdell 
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n'(134) = 45% ± 5% (peak) 

K(134) £ 2 radians. 

The effieicney at several different frequencies has been calculated from the 
data set. The results are given in Table 1. 


f (GHz) 
e 

ri(f e )%d; peak 

15 

67. 15 i . 13 

35 

63.73 * .67 

70 

60. 1 *2.3 

94 

55. 0 *3.6 

100 

53. 6 ± 4. 0 


Table 1 


Theorem 1 works very well so long as the peak-to-peak phase error 
does not get too large. For peak-to-peak phase errors above 3 radians at 
the measurement frequency the mathematical error of approximation 
increases so rapidly that the method becomes useless. It is interesting 
that this breakdown occurs at about 3/5 of the gain maximum frequency for' 
a reflector governed by the Ruze model. The gain maximum occurs when 
the rms phase error is about 1 radian or when’the peak-to-peak error is 
about 5 radians. Thus the theorem applies throughout what is considered 
the useful bandwidth of the antenna. On the other hand, the bound on the 
mathematical error is rather conservative since it bounds a weighted 
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average by its peak-to-peak value. One could actually perform the integral 
that this term bounds if the function L(x) were available from a STAIR pro- 
gram or a field measurement. Extreme accuracy would not be required 
since this integral represents only a small correction. If the appropriate 
integral were performed the method could probably be used with an efficiency 
measurement above the frequency of the gain maximum. 

Theorem 1 may be extended to the Cassegrain antenna case by con- 
sidering only the most significant terms. A rigorous result seems impossible 
because of amplitude and phase of the electric field on the primary reflector 

will inevitably be functions of frequency. The usual practice is to scale the 

f73 1 

primary feed horn and leave the subreflector unchanged . We can 
modify Theorem 1 to apply to a Cassegrain antenna by the following proce- 
dure. 

Let us define the per unit change in electric field from the measure- 
ment frequency to the estimation frequency as e{x) i. e. , 

E(f , x) =: E{f ,x) [1 + e(x)]. 

Also let us define the phase error produced by the feed system as .{.(f, x). 

It is important to distinguish this type of phase error from the phase errors 
caused by reflector surface deviations, 6(x), because -t(f, x) will not scale 
with frequency. Furthermore, suppose that the design efficiency which 
includes L-type phase errors is defined as ri o (f). 

With these definitions bounds similar to the ones given in Theorem 1 
are derived in Appendix C. This is done under the assumption that terms 
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up to third- order in e(x), -l(x), and 6(x) are retained. The only fourth order 

4 

terms retained are § (x) terms since the surface tolerance, 6{x), is small 
but more significant than e(x) or -L(x). In addition several terms were 
evaluated under the Ruze statistical model *. These bounds on the efficiency 
are given by 

’KV s W - PC V f 1 “ Gr) 2 +T io (f m )sup s+PB (r ) 2 

m m 

+ [Ti (f ) - n (f )] 0 2 (f ), for all P 
o m o e e 

and 

Tlty i VV - p t%y • Of ) Z + \ taf S 

m 

- PB ( tJ + C 1o (£ m > - "o< £ e» n2 < f e>- a11 p 

m 

where a(f^) the rms phase error in the aperture plane. 

As an example, suppose one wanted to estimate the efficiency of the 
JPL 85-foot antenna at S-band (2. 388 GHz) using X-band (8. 448 GHz) data. 
The necessary information is given by 

o(S) = 0. 166 radians 

K(X) = 2.9 radians'* 75 ^ 


*The statistical model is used to make the problem tractable. However, it 
is felt that the terms on which the statistical analysis is used will always be 
small. In fact, the only non-zero term is C T lo( ; ^m)“ r !o(^e)]0 2 ( : fe)* which is 
usually quite negligible. Thus, the method is only weakly dependent on the 
statistical model. 
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T| o (S) . =63.0% t?6 ^ 

ri (X) = 64. 7%^ 77 ^ 

T|*(X) = 44.7 ± 4. 9%{3 a )^ 78 ^ 

The efficiency estimate from the above bounds is 60.4 ± 1. 4%(3a). This 

f79 1 

compares favorably with the measured value of 59. 8 ± 8.9%(3cr) 

This data assumes that the error in the measurement is a random 
variable. The error bars are given as a percentage. Then the number of 
sigma units that this percentage represents is given in parenthesis. This 
procedure is probably the one most commonly used for efficiency measure- 
ments. The error is usually assumed to be Gaussian so the number of 
sigma units represents the confidence level for the measurement. For 
example, a Zp measurement means that the probability is 95% that the 
true value lies within the stated error bars. 

The error bars, stated in this way, may be used as peak errors as 
required by Theorem 1 and Theorem Z. Theorem Z will be presented in 
the following section. Then, the probability that the true efficiency lies 
between the inferred bounds is the greater than or equal to the probability 
that the measured efficiency lies between its error bars. 

Theorem Z is a generalization of Theorem 1. It gives a method of 
predicting the efficiency of a prime focus reflector antenna at frequency f 
based on measurements at two other frequencies. Here again one source 
of error is the error in the original measurements, but the mathematical 
error is considerably less since three terms of the cosine series are used. 
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Thus, it is possible to use a measurement well beyond the gain maximum 
frequency and obtain more leverage in reducing the measurement errors 
below gain maximum frequency. 

Theorem 2 

If the feed system of a prime focus reflector antenna is scaled at 
frequencies f , f , and f , and, 

G 1 M 

h'^) - rityi = Bj, 

h'(f 2 ) - r,(f 2 )S ^B 2 , 

j6(f e> x) - §(f e j y ) ! ^ K (fJ all x and y. 


and A 1 ^ and A^ are roots of the simultaneous equations 

’V, + A i' f l 2 + A 4 J 1 4 - ^l, ■ 


tv, + A 2 f 2 2 + A i f , 4 = 


4 2 


Then, 


Tl( f e ) <; n Q + PAJ f^ + QA^f g 4 + |F 1 |B 1 + |F 2 [B 2 + p o sup S, all P and Q 


and 

’l< f e )2T lo +J?A 2 f e 2 + QA 4 £ e 4 ' 1 F 1 1 B 1 ~ l F 2 I E 2 + H, taf S - 211 p =“ d Q " 

where 

f f 

S = |z:z = F 2 T[~ u]-F i T[j-u] + T[u]-( 1 -P) u Z /2! + (l-Q)u 4 /4‘; 
e e 


0 <; u £ K(f e )} 


and 
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£ e < P l Z ' Q£ e 2 > 


F„ = 


£ >‘, 2 -«. z > 

£ 2 2 * £ 2 2 - £2 > 


The remainder function T£x] is 


T[x] = cos (x) - £1 - x^/2i + x 4 /4I ]. 


This Theorem is useful when the peak-to-peak phasing error is 
large. Consider the case of The University of Texas 16-foot antenna as 
delivered by the contractor. It had an astigmatic phasing error as dis- 
cussed in Chapter III. This will be discussed in Chapter VI. Suppose one 
wanted to predict the efficiency at 35 GHz based on measurement at 70 


[80] 


GHz and 94 GHz. The data is repeated as follows, Cogdell 

p =65.0 

'o 

p'(94) = 44. 0 * 7 % (3 a ) 

p'(70) = 58. 0 ± 7% (3 a ) 

K(94) =4.3 

With this data the predicted efficiency at 35 GHz is 61. 58 ± 3. 11% (3a). The 
measured value is 62 ± 9%(3a)^^« This measurement is based on an 
astronomical measurement and assumes a disk temperature of 151° for 
Jupiter at 35 GHz. 
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The advantage of Theorem Z over Theorem 1 is that the mathematical 
error of approximation is considerably less for a given peak-to-peak phase 
error. Thus, the method is appropriate at higher frequencies where there 
is significant tolerance loss. The disadvantage of this method is that two 
measurements of the efficiency must be made. The method does, however, 
modestly improve the accuracy of the inferred efficiency. 

If the measurement errors have the same sign, the errors in the 
inferred efficiency tend to cancel. This fact can be exploited to obtain 
even more accurate inferred efficiencies. Whether or not there will be 
common errors in both measurements depends, of course, on how the 
efficiency is measured. However, if the same equipment and calibration 
techniques are used at both measurement frequencies some errors should 
appear in common. 

For example, in the efficiency measurement method used by Bathker^^ 
the uncertainty in the directivity of the standard gain horn is the dominant 
error. The directivity of this horn was calculated by integrating the experi- 
mentally measured pattern. If this procedure were repeated at two fre- 
quencies, the error in the horn calibration would tend to have the same 
sign. If the standard gain horns for the two measurements were scale models 
of each other, the error in horn calibration should be the same in both cases. 

One possible way of dealing with the common error is to assume that 
the measurement errors are correlated random variables. This assump- 
tion will be avoided, however, since these measurement errors are not 
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random variables, as pointed out by Cogdell^ and Ludwig^^. If they 
were true random variables, repeating the measurement should increase 
the accuracy. However, it is known from experience and- intuition that 
this averaging effect does not occur. Instead, errors tend to be systematic 
and bias the results. 

Instead of assuming that the measurement errors are correlated 
random variables, let us define the common error to be' c. 

Thus, 

■n*^) • ri( f x ) = c + v 1 
ri*{f 2 ) - r|(f 2 ) = c + r 2 

where r and r are the remainder errors which are not common to both 
measurements. Bounds for the efficiency can now be derived in terms of 
the peak values of c, r^, and r^, which will be denoted by C, R^, and R^. 

With this notation the bounds in Theorem 2 become 



all P and Q, 


and 



all P and Q, 

where all of the symbols are the same as before and 
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Suppose that 5% of the 7% error in the example considered before was 
in common. It should be emphasized that this is a hypothetical example. 


% 

= 65.0% 

r](94) 

= 44. 0% 

n(70) 

= 58. 0% 

R 1 

= 2% 

R 2 

= 2% 

c 

= 5% 

K (94) 

= 4. 3 radians. 


The bounds at several frequencies were optimized on a digital computer. 

The results are shown in Table 2. 

This table has two interesting features. The first is that the bounds 
at the measurement frequencies, 70 and 94 GHz, are much better than the 
original measurements. This happens because more information is used 
than is contained in the single measurement. The second interesting 
feature of this table is that Theorem 2 predicts a negative upper bound at 
200 GHz. The efficiency of an antenna is always positive, so a negative 
upper bound is a contradiction. It simply means that the data is inconsis- 
tent. In particular, this inconsistency can be traced to the fact that the 
measured efficiency at 94 GHz is too far below the measured value at 70 GHz. 
Thus, the parameter A is negative while it should be positive. 

This contradiction is not surprising since there is no reason to 
believe that these two measurements had any errors in common. The 
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Frequency Lower Bound Upper Bound 


GHz 

of Efficiency 
% 

P 

Q of Efficiency 

% 

P 

Q 

35 

60. 17 

1. 32 

1. 36 

61.61 

1. 23 

14. 02 

70 

51. 00 

1.00 

1. 00 

51.79 

1. 10 

3. 10 

94 

40. 50 

0. 94 

. 52 

41.73 

1.09 

1.70 

134 

18.63 

0. 80 

.22 

21.76 

0.99 

0.94 

200 

-22.03 

0.58 

.07 

- 5.98 

0.93 

0.73 


Table 2 

dominant error in these measurement is the uncertainty in the gain of a 
standard gain horn. The gain for the horn was measured at 94 GHz while a 
calculated value was used for the 70 GHz efficiency measurement. This 
example makes it clear that Theorem 2 can be used to check the consistency 
of antenna efficiency data. 

E. The Design Efficiency r| o 

In the above it was assumed that the design efficiency, n , is not a func- 

'o 

tion of frequency. In practice the design efficiency includes ohmic losses, 
which are proportional to the square root of the frequency. However, these 
effects usually are small. For example, at 100 GHz an aluminum reflector 
absorbs only 0.1% of the incident power. Thus, the assumed constancy of 
r) o with frequency will introduce negligible errors in the results. 
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Clearly the accurate calculation of r) is crucial to this method of 
scaling. The uncertainty in n was not included since it is typically much 
less than the uncertainty in the efficiency measurement. For a prime- 
focus instrument with solid supporting spars -q may be calculated within 
± 0.2%, see Ludwig^^* [87]. ^ would be a simple matter however to 

account for uncertainty in r| in the error bounds. 

F. Conclusion 

Two methods for predicting the efficiency at frequencies different 
from the one at which it was measured have been presented. The first 
method is based on a single efficiency measurement; the second uses two 
measurements. It has been shown that the error made in predicting the 
gain at a lower frequency is often considerably less than the measurement 
error. In fact is possible to achieve much better gain measurements by 
this method than any other known to the author. In addition, the theorems 
offer an easy way to check the consistency of antenna efficiency data. 
Unlike many antenna tolerance results, the method is independent of any 
statistical modeling in the case of prime-focus instruments, and is only 
weakly dependent upon it in the case of a Cassegrain antenna. 



Chapter VI 

Calibration Program for the 16- Foot Antenna 
In 1963 The University of Texas acquired a 16-foot reflector antenna 
designed to operate up to 14d GHz. It was first located in Austin, Texas, 
but in 1967 it was moved to Mount Locke, near Fort Davis, Texas. In 
Austin the antenna worked well at 95 GHz; however, after moving the 
antenna performance was unsatisfactory at this frequency. At 95 GHz 
no satisfactory focus position could be found, and the beamwidths were 
always broad in at least one of the principle planes. In addition, one side- 
lobe appeared at a level of -7 dB when the antenna was adjusted for maxi- 
mum gain. 

This chapter is an account of how the performance of this antenna 
was improved and a report of the subsequent careful evaluation. The 
antenna evaluation work reported here provides concrete examples of 
the methods presented in Chapters IV and V. In addition, this chapter 

rgg 

is meant to be a sequel to an earlier antenna evaluation report by Cogdell L 
The section A contains a detailed account of the measurement and 
correction of astigmatism in the reflector. The section B is a discussion 
of the patterns of the antenna. The patterns before the reflector adjust- 
ment are compared with the patterns after the adjustment. Both sets of 
patterns are compared with the theoretical patterns. Section C .is a dis- 
cussion of a careful antenna efficiency measurement at 134 GHz. A table 
is also presented which shows the inferred efficiencies at other frequencies 
of interest. Section D is a discussion of two antenna stability tests. 
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A. Adjustment of the Reflector 

The 16-foot antenna was known to exhibit astigmatism because it 
satisfied all of the properties given in Chapter IV. Thus, in June of 
1969 it was decided to adjust the figure of the reflector by shimming the 
backup structure. 

The positions of the shims that were used are adjacent to the eight 
bolts shown in Figure 25. The shim positions were numbered as Wl 
through W4 and El through E4 as shown in Figure 25. This figure also 
shows a top view of the antenna and the reference coordinates in the aperture 
plane. The coordinates x^ and u^ are in the declination plane while x^ and 
u are in the polar plane. Driving the antenna south increases the angular 
coordinate u 1 while driving the antenna east increases u . 

The adjustment was carried out by a two step procedure. The first 
step was making the cp direction agree with the servo coordinates; the 
second was reducing the magnitude of the astigmatism. The astigmatism 
was measured by observing the cp direction beamwidths. Thus, the cp di- 
rection was first made to correspond with the servo coordinates so that 
the cp dir ection beamwidths could be observed more easily. 

In each of the two steps outlined above the quantity of interest was 
first measured. Then, an arbitrary adjustment was made and the quantity 
of interest remeasured in order to get a feeling for the scale of the change 
produced by the added shims. This iterative procedure was continued 
until satisfactory results were obtained. For the cp direction adjustment 
of two different shim positions excluding the original state were tried. 




Fig. 25 
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The magnitude of the astigmatism was luckily reduced to an insignificantly 
small value by the first shim change. 

The initial cp-direction was known because it was calculated care- 
fully through the computer minimization procedure discussed as an 
example in Chapter III. However, it was decided to plot a contour map 
to see whether any difficulty was encountered in inferring the symmetry 
direction directly from a contour map. This first contour map is shown 
in Figure 26. The symmetry direction that was inferred from this was 
cp = -18°. This compares favorably with the value of cp = -9. 0° which was 
obtained from the computer minimization. 

In order to bring the major -axis of the elipse closer to the polar 
direction a shim of 0'.'020 was put in positions W1 and E3. Subsequently, 
another contour map was taken at approximately 10 dB down from the 
beam peak. The results of this are shown in Figure 27. The major 
axis of the best fitting elipse, drawn by eye, was taken to be +31°. It 
is apparent that the error was overcorrected with the 0'.'020 shim. 

Having removed the 0'.'020 shim a 0'.'008 shim was put under W1 and 

t 

E3. This value was chosen by linear interpolation between the -18° and 
+ 31° values. The resulting contour map is shown in Figure 28. From 
Figure 28 it is apparent .that the cp direction agrees with the servo co- 
ordinate direction, so the first step was complete. 

The magnitude of the astigmatism was then measured. The polar 
and declination beamwidths were taken at five different feed positions. 




POLAR ( Mi Hi - Degrees) 



Fig. 27 
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The results are shown in Figure 29. From this figure it is apparent 
that most of the astigmatic phase error still remains as might be ex- 
pected. The error was estimated as or = 1. 3. 

The action taken was to insert a shim of O'.'OIO under W4, Wl, E2, 
and E3. The astigmatism was then remeasured by taking beamwidths at 
four different shim positions. Since the beamwidths did not change signi- 
ficantly it was concluded that fortunately the astigmatic phase error was 
substantially corrected. Figure 30 shows a contour map taken after this 
adjustment was completed. It shows no significant directional preference. 

B . Patterns of the Antenna 

The method used for taking the patterns of the antenna has been 
discussed in detail by Cogdell^^. The antenna servo dri ves the antenna 
through the direction of the transmitter sight. At the same time the re- 
cording chart is marked every 0. 010 degrees by a trigger placed in the 
servo display network. 

Figure 31 shows a contour map of the error free pattern. The horn 
feed has the dimensions that correspond to the maximum forward gain for 
an antenna of F/D = 0. 5. These dimensions 1.25 X x 0.9 X. The efficiency 
achieved by the feed design is 65. 02%. 

Figures 32 and 33 show the antenna patterns reported by Cogdell. 
These were taken prior to the reflector surface adjustment discussed in the 
previous section. These figures show the design pattern or phase error 
free pattern drawn on the same scale. Figure 29 shows the 35 GHz pattern. 
It is very close to the theoretical pattern although one can see the beginnings 
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of pattern degeneration. Figure 31 shows the 95 GHz pattern. It is 
clear that from this figure that the observed pattern is much broader 
than it should be. The first sidelobes are greatly enhanced while the 
sidelobes after the first are relatively lower. As discussed in Chapter 
III this sidelobe arrangement is evidence that the phase errors decorrelate 
over distances comparable to the diameter of the antenna. 

Figures 34, 35, and 36 show the antenna patterns after the adjust- 
ment discussed in the previous section. Figure 34 is the observed pattern 
at 35 GHz plotted in the same scale with the predicted pattern. These two 
patterns agree very well. This means that there is very little phase 
error at 35 GHz. In addition, the excellent agreement of the observed 
and predicted pattern in the sidelobe regions at a frequency where there 
is little phase error means that the feed theory describes the illumination 
function very well. Figure 35 is the observed and predicted patterns at 
94 GHz. In comparison to Figure 33 this pattern shows that the phase 
error was much reduced by the adjustment. 

Figure 36 is the observed and predicted pattern at 134 GHz. At 134 
GHz the performance has deteriorated but the antenna is still very good. 
For some applications the overall efficiency of the radio telescope is not 
as important as the main beam efficiency. Since the beamwidths are 
broader than the design beamwidths by about ten percent in each plane, 
the energy that is lost from the peak direction is scattered near the main 
lobe. Consequently, the main beam efficiency is still quite high at this 


frequency. 
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C. Efficiency Measurement 

The efficiency of the antenna has been measured at 134 GHz, 
which is near the highest frequency of operation at the present time. 

As was shown in Chapter V this one measurement is sufficient to define 
the efficiency at all lower frequencies. Since the error bars scale 
down with the frequency squared, the inferred efficiencies at 15 and 
35 GHz have extremely close tolerances. This opens the possibility 
of making extremely accurate radio astronomical measurements at 
these frequencies. 

1. Measurement Technique 

The efficiency measurement was carried out by comparing 
the gain of the 16-foot antenna to a conical standard gain horn. The 
comparison was accomplished by disconnecting the feed from the 
radiometer input flange and connecting the standard gain horn to it. 

Thus, the measured gain is relative to the receiver input flange. 

The circuit that was used for the comparison is shown in 
Figure 37. One begins by finding the peak radiation level of the conical 
horn by taking the pattern of the horn in two orthogonal planes. Then 
with the horn pointed in the peak direction, a level was established on 
the recorder. After a few minutes, the receiver was switched to the 
16- foot antenna. Sufficient attenuation was then inserted on the precision 
IF attenuator to make the levels agree. The reading of the IF attenuator 
is the difference in gain between the horn and the 16-foot antenna. Sub- 
sequently the standard gain horn was reconnected to check for drift. 




GAIN COMPARISON CIRCUIT 


Fig. 37 
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This circuit and procedure is believed by the author to be 
the best available for the gain comparison measurements. The principle 
advantage of this method is that the gain difference between the standard 
gain horn and the large antenna is determined very accurately through 
the precision IF attenuator. The only disadvantage to this circuit is 
that it assumes linearity of the mixer over a 40 dB range. However, 
the mixer has been checked experimentally for saturation with negative 
results. The saturation effects begin at least 20 dB higher than the 
highest level encountered in the gain comparison measurements. 

The only reasonable alternative to this circuit is to use RF 
attenuation at the transmitter sight. This method also requires some 
mixer linearity because the RF attenuator must be calibrated by the IF 
substitution method. In addition, it has proved impossible to repeat 
calibrations of the RF attenuator at 134 GHz. This is consistent with 
a history of unreasonable efficiency measurements given by the RF 
attenuator method. Thus, the RF attenuation method was discarded 
as unuseable. 

The calibration of the standard gain horn is another crucial 

consideration in the overall accuracy of the measurement. The value 

used for the directivity is a value calculated by a method given by 

[gi] 

Cogdell . This method gives values which are in excellent agreement 

with several horns of varying dimensions measured experimentally by 
[92] 

King. Some had considerably larger flare angles than standard gain 

horn used for this measurement. The major shortcomings of Cogdell's 
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method is in accounting for the effects of the non- zero horn flare angle 
Thus, it is felt that the calculated value for the horn gain is very close to 
the actual value since Cogdell's theory seems to work well under con- 
siderably more demanding circumstances. 

It is difficult to put error bars on a measurement of this kind 
because all errors tend to be systematic, so repreating the measurement 
does not increase accuracy. The procedure that has been adopted is to 
give the accuracy of each step in the measurement. Then the total error 
is computed by assuming the errors are random variables and the error 


limits given are values of 2a(e) where e is the error. This procedure is 

us til 
[95] 


[94] 

not really justifiable as pointed out by Ludwig , but it is the customary 


procedure. 

2. Results 

The results of the efficiency measurements are given in 
Table 3. All of the values are given in dB. The 16-foot antenna gain is 
simply the gain of the horn plus the difference in gain registered on the 
IF attenuator. The gain of the standard horn at the receiver input flange 
is the directivity of the horn minus the omic loss in the horn and the loss 
in the 2 foot waveguide connecting section. 

The efficiency is the gain of the 1 6-foot antenna minus the 


gain for 100% efficiency. The tolerance loss which is the loss from phase 
errors is the directivity of the antenna minus the design directivity. The 
directivity is higher than the design gain by the omic loss in the feed. 
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Table 4 gives the efficiencies at other frequencies calculated 
from the Cassegrain antenna version of Theorem 1 given in Chapter V. 

It is necessary to use this extended version of Theorem 1 because the 
actual feeds which are used are not perfectly scaled. The theoretical 
efficiencies were calculated with the program given in Appendix A using 
the feed dimensions given in the table. The feed omic loss factor be- 
comes appreciable at higher frequencies. Since this loss does not scale 
with frequency, it was taken out before the scaling was done. Subsequently, 
the feed omic loss was reinserted to calculate the total efficiency which is 
relative to the receiver input flange. 

D. Antenna Stability 

The observed patterns and gain measurement in the previous sections 
were taken in the hours from 2 a.m. until sunrise. This time period has 
been found to be best because the atmospheric induced signal fluctuations 
are smallest during these hours. Since these measurements were done 
under more or less ideal thermal conditions, it is mandatory that the 
antenna stability be tested under less favorable conditions. 

There are two main effects that might occur. One of these is 
differential gravity loading. The antenna is always pointed in the direction 
of the transmitter site to make gain and pattern measurements. Thus as 
pointed out before gravity loading might cause astigmatism in the antenna 
when it is pointed in another direction. 

Temperature effects are the other possible source of instability. 
Strong temperature gradients, especially, might affect the antenna 
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Summary of Gain Measurement of the 
1 6- Foot Antenna at 134 GHz 


Item 

Value (dB) 

2a (dB) 

A. 

Gain Difference 

40.55 

0. 10 

B. 

Gain of Conical Horn 

32.55 

0.27 


• "f" 

1. Directivity of Horn^ 

34. 98 

0. 25 


2. Omic Loss in Horn 

0. 1 

0. 1 


3. Loss in Waveguide 

2. 33 

0. 05 

C. 

Gain of 16- Foot Antenna 
at 1 34 GHz 

73. 10 

0.28 

D. 

Gain for 100% Efficiency 
at 1 34 GHz 

76. 71 

— 

' E. 

Efficiency of Antenna 

-3. 61 

6. 27 


at 1 34 GHz 

(43.6%) 

(±2. 8%) 

F. 

Design Gain 

74.55 

0. 1 


1. Design Directivity (65. 02%) 

74. 85 

-- -- 


2. Feed Omic Loss 

0. 3 

0. 1 

G. 

Tolerance Loss 

-1. 45 

0.29 


Calculated by Cogdell's method as discussed in text. 

•A- 

Estimated 

Measured radiometrically 


Table 3 



EFFICIENCY CHARACTERISTICS OF THE U.T. 16 FOOT ANTENNA 

JANUARY 1969 


FREQUENCY 

FEED DIMENSIONS 

FEED 

design 

EFFICIENCY 

WITHOUT 

FEED 

LOSS 

(%) 

TOTAL 

TOTAL 

TOTAL 
EFFECTIVE 
AREA 
( m 2 ) 

(GHz ) 


a '\ 

b 

(inches) 

b /\ 

LOSS 

(dB) 

EFFICIENCY 

(%) 

EFFICIENCY 

(%) 

GAIN 

(dB) 

15.0° 

0.986 

1.25 

0.709 

0. 900 

0.01 

65.02 

64.73*0.0# 

64.58*0.09 

55.79 ± 0.01 

12.061 0.02 

35.0° 

0.422 

1.25 

0.304 


0.03 

65.02 

63.50±0.46 X 

63.06 ± 0.46 

63.04*0.03 

11.781 0.09 

35.0 R 

0.400 

1.18 

0.285 

0.844 

0.03 

64.57 

63.05* 0.4# 

62.61 ±0.46 

63.0I±0.03 

11.70*0.09 

70.0° 

0.2H 

1.25 

0.152 

0.900 

0. 10 

65.02 

59.2 ± 1.7 X 

57.8 ± 1.6 

68. 69 ±0.11 

10.8010.30 

94.0 R 

1.164 

1.30 

0.113 

0.899 

DB 

64.92 

55.0 ± 2. 5 X 

53.1 ± 2.4 

70.88±0.I9 

9.92 ± 0.44 

I00.0 R 

0.164 

1.39 

0. 113 

0.957 

0. 17 

64.06 

53.0 ± 2. 6 X 

50.9 ± 25 

71.23 ±0.2i 

9.51 1 0.47 

134. 0 R 

0.104 

1.18 

0.076 

0. 862 

0.30 

64.68 

46.6 ±3.1 

43.6 ± 2.9 


8.13 1 0.53 

I40.0 0,R 

0.104 

1.24 

0.076 

0.901 

0. 30 

65.02 

45.7 ± 3.6 X 

42.6 ± 3.3 

— 
73.38 + 033 

7.96 1 0.62 


All confidence levels are 2 cr levels in ihe same units as the quantity . 
0 Feed designed for maximum forward gain 

D 

Feed on Rathlon radiometer 
* Measured 

x Based on Theurem t, ChapterTT 


TABLE 4. 


121 





















































































122 


structure. A strong temperature gradient is produced, for example, 
by the sun shining on part of the antenna and the other part is in shadow. 

The first experiment performed was an experiment to check for 
temperature gradient effects. Just before sunrise the antenna was 
pointed at the transmitter with the chart recorder calibrated on a rather 
sensitive scale. Then the system was left untouched for about an hour. 
During this time the sun shone on more of the antenna. The experiment 
was terminated when the sun shone on about half the antenna. The record 
showed no evidence of gain deterioration. 

This experiment is significant because it exposes the antenna to the 
worst possible conditions of differential heating. The first effects of 
differential thermal loading of the structure is to steer the beam away 
from the equilibrium position. Since no significant gain reduction occurred 
there is apparently no significant thermal loading. 

As mentioned earlier differential gravity loading should produce 
astigmatism. It was known that there was no astigmatism when the 
antenna was pointed toward the transmitter sight, but there could well be 
astigmatism at some other pointing. Thus, the second experiment was 
designed to detect astigmatism at other pointings. This was done by the 
solar limb crossing experiment discussed in Chapter IV. From this data 
one can deduce the beamwidth of the antenna in the direction the antenna 
crosses the limb. Then as discussed in Chapter IV the astigmatism can 
be determined from plotting the beamwidth vs focus position. 
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For this experiment, thirty- six polar and thirty- six declination 
scans were made of the sun. Figure 38 shows a typical scan. The 
experimental procedure was to make two polar and two declination 
scans. Then, the axial feed position was changed and another set of 
scans were taken. Three feed positions were used, with twelve scans 
in each direction for each feed position. 

Table 5 shows the results of this study. This table shows no 
significant astigmatism. In addition, the beamwidths at the shim 
position nearest the optimum agree very well with the value measured 
on the pattern range. This occurs even at an hour angle of +42° and a 
declination of -20°. This is a very different pointing from the trans- 
mitter site which is at an hour-angle of -69° and a declination of -34° 
From this data it was concluded that there was no significant gravity 
loading of the antenna structure. 

E. Summary 

This chapter illustrates the theory that has been developed in 
Chapters IV and V. It is a detailed account of how this theory was put 
into practice. The work was successful in substantially improving the 
performance of the 16-foot antenna and establishing its stability in all 
pointing directions. In addition, the efficiency of the antenna has been 
accurately measured at 134 GHz, and the efficiencies at other frequencies 
of interest inferred by the method given in Chapter V. 



RELATIVE SCALE 
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Declination and Hour Angle Beam widths in Milli- Degrees 


Hour Angle in Degrees (x 10 ) 



335 

12 

42 

Pattern 

HA 

33.2 

33. 0 

32. 3 

33. 5 

DEC 

34. 5 

33. 3 

34. 5 

32. 5 


T/X = 

0.23 




Hour Angle in 

-3 

Degrees (x 10 ) 


312 

354 

21 

Pattern 

HA 

3-7. 6 

36. 0 

36. 0 

— 

DEC 

38. 4 

35. 0 

34. 6 

— 


v/x = 

0. 68 




_ 3 

Hour Angle in Degrees (x 10 ) 


320 

2 

31 

Pattern 


HA 

— 

46. 8 

45.0 

— 

DEC 

40.2 

43. 9 

46. 0 

— 


T/X = 

1. 13 




Declination of sun -20° 

All values except the pattern range measurements are ±1. 5 milli-degrees 

(2a) 


Table 5 



Appendix A 


An Antenna Radiation Pattern Program 
Using the Fast Fourier Transform 

A. Introduction 

Calculating the far field radiation pattern of a reflector antenna 

from the aperture fields is a problem which has attracted considerable 

interest. With ordinary analytical methods one is restricted to a very 

small class of aperture illuminations and phase functions. For example, 

one can calculate the radiation pattern of a circular uniform phase 
[ 97 ] 

aperture with illumination proportional to 


1 - 


r 

KTJ J 


-, P 


o 

where r^ is the radius of the aperture. This form is convenient because 
its Fourier transform is a Bessel function. This aperture distribution is 
also useful for tutorial studies of the general effects of tapering on reflector 
antennas 

The engineer involved in the evaluation of a large reflector antenna 

has a need to calculate the radiation pattern in a much more general con- 
[ 99 ] 

text. Spencer has given an infinite series expression (later corrected 
by Fagot f or the radiation pattern of an antenna with phase errors. 

However, this series has two disadvantages. The first is that it requires 
one to calculate the derivatives of the design pattern with respect to the 
angular coordinates. The second is that no convergence criterion is given 
for the series. In order to treat realistic problems Ludwig 
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[102] [103] 

Dion , Bao , and others have turned to numerical methods. By- 

calculating the radiation pattern numerically one is not restricted to any 

specific tapering or phase error function 
. [104] [105] 

As Dion and Bao have pointed out the computation required 

is enormous. Suppose one wanted the radiation pattern calculated at an 

2 

(NXN) grid of points. This requires the calculation of N two dimensional 

integrals. If each integral is approximated by an NXN sum then one must 

4 

carry out approximately N additions and multiplications. This is too 
large to be practical. For example, with N = 128 on the CDC 6600 this 
procedure would take over one hour. 

In this Appendix a computer program will be presented that takes 
advantage of the efficiency of the Fast Fourier Transform (FFT)^^^. 

The FFT can be used to perform the same calculation as outlined above 

by using the FFT 2N times. The FFT itself requires approximately 

* [107] 

N log N additions and multiplications . Thus, the time savings 

factor is 


r = 


2 N • N log N log N 

c* 


N' 


N 


.-3 


If N .= 128 then r < 10 . Thus the FFT method is more than a factor of 

one thousand more efficient than the direct method. This efficiency makes 
possible the model fitting discussed in Chapter III. This method also uses 
storage very efficiently because only a single N x N complex array is 
required. 
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A program similar to this one has been developed by Pratt and 
Andrews^" . Their program calculates the defraction pattern of 
light from an arbitrarily shaped aperture in a plane screen. However, 
it would be difficult to use this program to calculate an antenna pattern 
because of the blockage problem which will be discussed later. In addi- 
tion, the program given here allows more flexibility in choosing the grid 
size for both the aperture plane and the far zone angular coordinates. 

In the following section the mathematical development of the program 
is given. In Section C the program itself is discussed and a listing is 
given. In Section D the computed results for a uniformly illuminated 
aperture are compared with the exact results. This program will calcu- 
late the radiation pattern of a circular aperture excited by electric field 
with arbitrary amplitude and phase. It assumes linear polarization of 
the aperture fields. The program allows one to account for central cir- 
cular blockage and the blockage of two spars at arbitrary angles to the 
reference coordinate system. The spars are, however, assumed to pass 
through the origin of the reference coordinate system. This assumption 
and the assumption of a circular aperture are not required by the mathe- 
matical development, but they are reasonable assumptions and simplify 
the program coding. 

B . Mathematical Development 

As shown in Chapter I the far field gain function can be written as 
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g(u i- u 2 ) = vl I<u r u 2’| 2 

where 


(Al) 


T# . rf T , , jkL(x r x 2> jk(u i x i + V X 2 > „ 

i {u r u 2 ) J J T ^ x i * x 2 e e 

— CO 

(A2) 

The tapering function Tfc^.x^) is defined as the electric field in the aperture 
plane normalized by the square root of the input power and the free space 
impedance, i.e. 


E(x ,x ) 

T(x l’ x 2 i = /P. I 


m o 


jktfr^x ) 

e 


G f (9 7 , cp 7 ) 

Z 2T ~ 2 

^(f+^) 


jk't(x 1> x^) 
e 


Equation (A2) is particularly convenient for computer calculation 
because the feed characteristics can be thought of as an input and the far 
field gain pattern of the reflector as the output. In the case where the feed 
is an extended source, as it is in a Cassegrain feed, one must take explicit 
account of the phase shift, kt (x^ , x^), caused by the feed system. For a 
point source feed there is no phase shift caused by the feed, so -l may be 
set to zero. 

A typical tapering function is shown in Figure Al . These functions 
are gradually changing, well behaved functions, except where the antenna 
feed system blocks the aperture. In these blocked areas the functional 
value is zero. The blocked regions are characterized by having at least 
one dimension much smaller than the diameter of the antenna. 



TCX, ,X 2 ) 


tapering function 
WITH BLOCKAGE 
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This latter fact leads to a rather fundamental problem in choosing 
the integration grid size for the integral, I. If the grid size is chosen 
appropriate to an unblocked antenna, it will be much too course to account 
for the blocked areas effectively. On the other hand if the grid size is 
chosen small enough for the smallest blocked area, then the core storage 
and time required to calculate the integral is exorbitantly large. 

A solution to this problem is to break the tapering function into 
several parts so that it is written as the unblocked tapering function plus 
several corrections. Then the contribution of each component is calcu- 
lated separately and the results totaled. This method takes advantage of 
the linearity of the integral, and it may be thought of as the engineering 
method of "superposition. " The tapering function is written 

T (x) - T , <x)‘- T <x) - .T (x) - T (x) 

— ub — c — si— sZ — 

The function is the tapering function the antenna would have without 

blockage. The functions T , T , T are the blockage functions to account 

C S i S u 

for the central blockage and the blockage of two spars. These functions 

have a value of zero outside the area to which they refer and the value 

of T , inside this area. This process of decomposition is illustrated 
ub 

in Figure A1 . 

The components of the integral I that correspond to each of the 

tapering functions will be written I . , I , I and I _ after their 
c ub c si s2 

respective tapering functions. The technique used in the program is to 
calcuate the unblocked integral, I which is by far the dominant 
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contribution. The others are accumulated as corrections. This is a 
standard digital computer technique and conserves a great deal of 
storage. 

Calculation of the Unblocked Integral 

The three basic steps in calculating the unblocked integral, I 
are the calculation of the tapering function on an array of points, the 
introduction of a linear phase grading, and the performance of a row- 
wise and column-wise transformation of the array with the FFT. The 

result of this transformation is that I (u , u ) appears in the array at 

ub 1 2 

a grid of points. 'The tapering function is an input to the program and 
may be put in as a table of laboratory measurements. It may also be 
calculated from a theory for the feed system. The introduction of linear 
phase grading merely translates the points in the transformed array. 
This is done for convenience so the axial gain will be in the center of 
the array. The remainder of this section is devoted to the mathematical 
development of the FFT method of calculating this integral. 

The integral that must be calculated is 

00 

jkL( v * 2 > jk(u jXl + u 2 x 2 ) 

W'VVJJ T < x !' x 2> e 

CO 

The (x , x ) - coordinate system is the natural coordinate system for the 
reflector, but it proves inconvenient for Fortran programming. Thus', 
let us introduce a (y^Y^) “ coordinate system as 
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y l = *1 +h l 


y 2 * *2 + h 2 


The purpose of this translation is to shift all of the non-zero values of 
T into the first quadrant. Figure AZ shows these two coordinate systems. 
In (y , y ) - coordinates the integral above can be written 

1 u 1 




-jkfu^ 


u 2 h 2> 


U T(y l- 


T ’ y 2 


■v 


jkL(y L - 


h l' y 2 


h 2> 


jktujyj + u 2 y 2 ) 

e dy l dy 2' 

The computer calculation of an integral requires that' it be approxi- 
mated by a sum. A uniform grid of points that are Ay, and Ay? apart in 
the two coordinates are used, so that the integral may be approximated as 

+ u 2 h 2 ) * * 

1 ub (u r u 2 ) = e 4y i 4y 2 L L T(miy i • V nAy 2 ’ V 

m=l n=l 

jkLdmAyj - hj, mAy 2 - h.^) jk^mAy]^ + u 2 n AY 2 ) 
e e . (A4) 

In the following, the explicit form of T and L will not be written out, 

but will be written as T(m, n) and L(m, n). It will be understood that the 

indices (m, n) stand for the appropriate place on the aperture plane. Thus, 

Equation (A4) can be written as 
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WvV*® 


-jk{u 1 h 1 + u^h^) 


K K 

A y i&y 2 ^ ^ T(m,n) e 
m=l n=l 


jkL(m, n) 


jk(u 1 mAy 1 + Ta 2 n&y 2 ) 


(A5) 


All of the non- zero values of T(m, n) e jkL(m, n ) are pi ace( j an 
N x N complex array which is stored as a real array, PR, and an imag- 
inary array, PI, in PAT. N and K are unequal in general with N < K. It 
will usually prove convenient to choose 1C a factor of two to four larger 
than N, For the FFT routine used here N and K must also be powers of 
two. The reason for this will become apparent later. 

It is convenient to choose h^ and h 2 so that the value of T and L at 
the center of the aperture appear in the center of the array. This may 
be accomplished by choosing h^ = N/2 Ayj, and h 2 = N/2 Ay^* 

The u-variables in Equation A5 are continuous variables, so they 
must be rewritten as discrete indexed variables. Let u^ and u 2 be defined 


as 


u i = lu i 


u 2 = ‘ 


s - ±1 , 2, 
t - ±1,2, 


(A6)' 


where 


Au, = 


2tt 


1 K k Ay^ 


and 
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2tt 

K k Ay 2 ‘ 


The grid size for the u-array is chosen so that the integration grid 
size and the wavenumber fall out of the second phase term in Equation 
(A5). This normalized form is required by the FFT. Substituting the 
values discussed for h , h , u , and u in Equation A5, one has 

1^1 u 


.N , , . K K .2tt, , 

-3-.tt(s + t) _ _ xt , » Jyr ( ms + nt ) 

I ^(s,t) = e Ay l Ay zZ X T(m,n)e J ’ e 


ub 


m=l n=l 


(A7) 

the indices s and t represent a point in the angular coordinates u^ and 
u just as m and n represent a point on the aperture plane. 

Half of the s and t indices are negative, and since negative 
subscripts are not allowed in Fortran programming, s and t must 
be biased so that the only positive indices are used. This may be done 
by introducing a linear phase grading into the input which shifts the axial 
direction (s = t= 0) to the middle of the array. Let 


and 


/ 


s 


= s 


+ 


N 

2 


t 7 



With these translations Equation A7 becomes 
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■ N , , 

e 


+ t') 

A y i & y 2 


K K 

£ X T{m,n)e jkL(m ’ n) 


m^l n“l 


-jTT— (m + n) j-gr (ms ' + nt') 
e e 

Equation (A8) ma y be regrouped as follows: 


(A8) 


•N f , . A 


/ , - j K + 1 ) r V 1 - k z' r* 

I - (s ’ t) = e Ay l Ay 2lA e K.L 


ub' 


K .2n / K .2it . 

• ms nt 

K z' T K 

e 


P(m, nf)j 


m=l 


n=l 


(A9) 


where 


_ . , jkL(m.n) {m + n) 

P(m, n) = T(m,n) e J ' e 


The Fast Fourier Transform is a program that replaces an arra y , 
sa y B(n),with its discrete transform, i. e. , 


K 

B(n) -» ^ B(m) 


.2rr 

j^-m n 


m=l 


One can see that the calculations indicated in Equation (A9) can be done 
by first applying the FFT to each row of P(m, n) and replacing the results 
in that row. This procedure gives an arra y P(m, t'). Then the FFT is 
applied to each column of P(m, t') and the result replaced in that column. 
The result of this is the desired integral except for a phase shift factor. 
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Correction for Central Circular Blockage 

The central blockage correction can be written 

jkfu^ + u^x ) 

e e dx dx 

central 

region 

The region of the central blockage is small compared with the variations 
in T and L, so T and L will be assumed constant for this calculation. 


WVV “ 


J J 


T ub ( V*2> 


*cb^ U l ’ U 2> = T * 0, e 


jkL(0, 0) 


H 

central 

region 


jkOljX, + u 2 x 2 ) 


dx 1 dx 2 . 


Changing variables to r and § such that 


= r cos ^ 
x 2 = r sin ^ 

and 

u^ = 0 cos cp 
u^ = 0 sin cp, 

one has 


I cb lu l’ u l = T( °' 0) 


e jkL(0,0) 


JJ e jkSrco s (9- 5) rdrd5 , 

0 0 


where R, is the radius of the blocked area. The integral in 5 is just the 
b 

zero order Bessel function, so one can write 
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1 

I , 11 ) = T (0, 0 e j kL(0, 0} 2tt R* f J {k0 R z) z dz. 
cb 1 L b J o b 

0 

This last integral is also a Bessel function, i. e. 




(A10) 


where 9 may be calculated as the Pathagorian sum of u^ and u. 


2 

‘2 • 


9 “A i 2 + u ; 

When written in terms of the indices and t' , Equation (AiQ)becomes 


I cb (s',t') = T(0, 0) e ikL<0 ' 01 


where 


(All) 


G^)V - D" <^) 0' -D' • 

Spar Blockage 

The spar blockage correction may be written, using vector notation, as 


R, . 2 


jkL(x ,x ) T 

I S (U 1’ U 2 ) = JJ T ub (x l’ X 2 )e eJ “ “^1^2 (A12) 

spar 

Let us transform coordinates so that the primed coordinates agree with 
the direction of the spar as shown in Figure A2. The primed coordinates 
are given by 
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x' =Rx 


where 


R = 


cos 8 - sin 8 

sin 6 cos 6 


Also let us define the primed angular coordinates as 


n' =Ru 


jkL{x^ , x^) jkfu^Xj 
e e 

spar 


W 


u. 


;WJ 


T ub< X l' X 2> 


+ u„ 



dx^ / d X 2 ' ' 


(A13) 

The x ' integration is now trivial since I and L do not vary much over 
* ^ 

the small extent of x^*, Thus Equation (A13) may be written 


“ jkL(x ,x ) jku 'x ' 

I T ub ,x l' x 2> e « dx/ (A14) 



Now we need to change variables to y so that all of the non- zero values of 

T , are to the right of the origin. Let 
ub 


x ][ ' = y - h, 
so Equation (A14) becomes 

-i k u i' h Z j ku / y 

I s (U] '.u^') = w e J z(y) e dy 

- 05 

where z(y) is a complex quantity given by 
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jkL(x , x ) 

a(y) = T ub <X l ,X 2 )e 


where 


x^ = (y - h) cos 6 
x^ = (y - h) sin 6 . 


Approximating the integral by a sum of k terms one has 


-jk u ' h jk u, ' mAy 

I g {Ui '.u^') ~ w Ay e ^ z(mAy) e (A15) 

m=l 


Let 


Uj ' = v Au^ ' v = ± 1 , 2, 3, .... (A16) 


where 


/ 2tt 

An, - — < 

1 K k Ay 

Substituting Equation (A16) into Equation (A15) and letting h = — Ay. one 

c* 

obtains 


K ,2tt 

“Jttv h7 mv 

I g (v) = w Ay e ) z(m) e , (A17) 

m=l 

where it is understood that I (v) is I {u / ,u_ / ) where u. / = v An, ' . In 

s s 1 2 1 1 

addition, z(m) stands for z{mAy). 

Some of the v's of interest are negative. This is inconvenient for 


Fortran programming so let 
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v'= v + K/2. 

Thus Equation A17 becomes 

z K K ,2 tt t 

-j TT (V - 7 ) r-, _ iTTrn 3TF “V 

I (v*) = w Ay e ) z{m) e e (A18) 

m=l 

The sum in Equation (A18) may be calculated, with a single appli- 
cation of the FFT. The array z(m) is calculated first. Second one 
modifies z(m) with the linear phase grading e * 5TTm , Third, z(m) is 
transformed with FFT, and finally the transformed array is multiplied 
by the required constant factors and linear phase grading. 

The value of v / that corresponds to s ' and t y will now be calcu- 
lated. The angles u and u are given by 

1 u 

2tt 

u. “ Tr~, s 

1 K k Ay 1 

_ 2tt 

u 2 ~ K k Ay 2 * 

The value of u^ ' is then 


2tt 


u, = 


K k Ayj 


s cos 6 - 


2n 


K k Ay- 


t sin 6. 


Substituting Equations (Al6) for u^ ‘ into the equation above and cancelling 
the common terms, one obtains 


v — s cos 5 - E 

Ay x Ay 2 


sin 6 . 
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Now substituting v' , s' , and t' for v, s,t, one obtains 

v' = (s' - H) cos 6 - (t / - H) sin 6 + K/2. 

^ 1 ^2 

A slight problem arises in that the values of v 7 needed for the 
blockage correction are not necessarily integers. However since this 
correction is a small one, it is felt that interpolating linearly between 
integer values is sufficiently accurate, 

C. The Program 

The first input to the program is a card that has the diameter of 
the antenna in feet in the first field and the wavelength in meters in the 
second. All fields are E 10. 3. The next card specifies the angle between 
the far field pattern grid points in degrees. If this card is left blank 
the program calculates its own increment to optimize accuracy. This 
latter mode of operation is recommended unless the grid interval must 
be specified, as when it is being compared to measured data. The third 
card specifies the radius of the central blockage circle in inches. The 
fourth card specifies the width of the two spars in the first field in inches. 
It also specifies the angle the two spars make with the {x^ , x^)- coordinates 
in the second and third fields. 

The subroutine MAG(X1,X?) calculates the magnitude of the tapering 
function. The tapering function is assumed to be real in this program. 

The subroutine PHASE (Xl,X2 t Z, NX, PH) calculates the phase 
function kL(x). The phase in radians is returned as PH. The Z vector 
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is a parameter vector. The number NX is the maximum number of 
parameters. 

The internal program parameters KEXP and NEXP determine the 
grid spacings, accuracy, and storage requirements of the program. 

The unblocked integral is carried out over an N X N array of points 
with spacing Ay^ and Ay^> where 

NEXP 
N = 2 . 

The grid spacing of the points in angle space is given by Equation (A6) as 


- 2tt 

A “l “ K k iy 2 

and 

* - 2rr 

AU 2 K k Ay 2 * 

where 


(A19) 


K = 


KEXP 

2 


If the grid spacing in angle, space are left unspecified, then Ay and Ay 9 
are calculated to just fill up the N X N array. If kv. and Au_ are specified, 

X La 

on the other hand, Ay^ and by^ are calculated from Equation (A19). Since 
K and N may be different, the grid intervals may be adjusted somewhat 
for the particular application. Increasing N makes the program more 
accurate while increasing K makes the grid in angle space finer. In- 
creasing N, however, increases the program storage requirements. 
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One detail of the coding differs from the mathematical development 
given in the previous sections. The FFT routine used here actually shifts 
one grid spacing as well as transforming. That is, if B(m) is a complex 
array being transformed, then 

* {m - 1) { kV 1) 

B(K) - ^ B (m) e 

m=l 

D. Verification 

The calculation of the unblocked integral was checked by setting 
the blockage corrections to zero and calculating the pattern of a uniformly 
illuminated constant phase aperture. This was done by substituting another 
subroutine for FEED. FEED normally calculates the gain of the horn feed. 
This subroutine calculates the gain of the point source feed required to 
give a uniform aperture illumination, Silver^ For this test N = 128, 

K = 512, D = 16 feet, \ = 3.2 mm, and f = 8 feet. The result was an 
efficiency of unity to more than five significant figures and a normalized 
pattern given in Table Al. 

The central blockage and spar blockage corrections were checked 
by setting the tapering function to unity everywhere and printing out these 
corrections separately. Then several points were checked by hand. 



Patterns (dB) 


angle 


(xicr 3 ) 

Calculated 

Exact 

0 

0. 00 

0. 00 

10 

-0.77 

-0. 77 

20 

-3. 23 

-3. 23 

30 

-8. 06 

-8.06 

40 

-18. 32 

-18. 31 

50 

-24. 05 

-23. 98 

60 

-17.64 

-17. 59 

70 

-19. 68 

-19. 61 

80 

-30.50 

-30. 27 

90 

-28. 56 

-28.72 

100 

-23. 83 

-23. 91 

no 

-26.45 

-26. 56 

120 

-30. 96 

-30.74 

130 

-27. 98 

-27.86 

140 

-31.87 

-31.81 

Table Al 


Verification of Program 



no non non 
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E. Listing 


PROGRAM PAT<INPUT, OUTPUT) 

DIMENSION PR AND PI TO NXN. 

DIMENSION PR(128, 128), PI(128, 128) 

DIMENSION A AND B TO K. 

DIMENSION A( 1 024), B( 1 024) 

DIMENSION D AND E TO N. 

DIMENSION D{128), E(128) 

DIMENSION Z(15) 

DIMENSION SIG(2) 

DIMENSION MM(2), IDEC(2) 

DIMENSION IX(5), JX{5), MX(5) 

DIMENSION XX AND YY TO 2K. - 

COMMON/ SCRFFT/ XX(2048)/FFF/ YY(2048) 

INTEGER H 

INTEGER HK 

REAL LAMDA 

REAL MAG 

RES1 = 0.5 

RES2 = -0. 5**4 

RE S3 = 0. 5**7 /3. 

RES4 = -0. 5**7 / 3. /2. /4. / 3. /2. 

10 = 3H 0 
15 = 3H -5 
110 = 3H-10 
115 = 3H-15 
120 = 3H-20 
125 = 3H-25 
130 = 3H-30 
135 = 3H-35 

JX(1) = 10HTICK MARKS 
JX(2) = 10H ARE 5PI*D 
JX(3) = 1 OH /LAMDA 
MM(1) = 10HCONTOUR MA 
MM (2) = 1HP 
13 DR = 4H 3DB 
I10DB = 4H10DB 
I20DB = 4H20DB 
I30DB = 4H30DB 
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IPOL = 5H UZ 

IP AT = 7HPATTERN 

IDEC(l) = 10H U1 

PIE = 3. 1415926535897932384 

NEXP = 7 

KEXP = 10 

N = 2**NEXP 

K = 2**KEXP 

H = N/2+1 

READ 28, DIAMF, LAMDA 
DIAM = DIAMF *0. 0254*12. 

SCALE = PIE *DIAM/ LAMDA 
WAVENR = 2. *PIE/ LAMDA 
READ 28, DDl, DD2 
C 

C IF DDl IS ZERO THE PROGRAM 

C CALCULATES ITS OWN DDl AND DD2. 

C 

IF {DDl. NE. 0. 0) GO TO 1 

DU1 = 2.*PIE/K/WAVENR/DIAM*(N-2) 

DU2 = DU1 

DDl = DU1*180. /PIE 
DD2 = DU2*1 80. /PIE 

1 CONTINUE 

DU1 = DD 1 *PIE / 1 80 . ' 

DU2 = DD2*PIE/180. 

READ 28, RB 

READ 28, W, SIG(l), SIG(Z) 

PRINT 29 

PRINT 30, DIAMF, LAMDA 
PRINT 31, NEXP, KEXP 
PRINT 32, DDl, DD2 
PRINT 33, DU1, DU2 
PRINT 34, RB 
RB = RB*0. 0254 
PRINT 35, SIG(l), SIG(Z) 

DO 2 I = 1, 2 

2 SIG(I) = SIG(I)*PIE/180. 

PRINT 36, W 

W = W*0. 0254 

DY1 = 2. *PIE /K / WAVENR / DU 1 
DY2 = 2. *PIE/K/ WAVENR/DU2 
EXT1 = DIAM/DY1 
EXT2 = DIAM/DYZ 
PRINT 37, EXT1 
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PRINT 3 8, EXT 2 
IF (EXTi .GR. (N- 1)) STOP 1 
IF (EXT2 .GT. (N-l)) STOP 2 
C 

C THIS COMPLETES THE INITILIZATION PROCEDURE. 

C 

C 

C CALCULATION OF THE UNBLOCKED INTEGRAL FOLLOWS. 

C 

CALL SECOND (QT) 

CONST = PIE*N/K 
DO 3 I = 1, N 
DO 3 J = 1, N 
XI = (I-H)*DY1 
X2 = (J-H)*DY2 
TMAG = MAG(X1, X2) 

CALL PHASE (XI, X2, Z, NX, PH) 

PH = PH- (I- H+ J- H)*C ONST 
PR(I, J) = TMAG*COS(PH) 

PI(I, J) = TMAG *SIN(PH) 

3 CONTINUE 

CALL SECOND (QT1) 

Q = QT1 - QT 
QT = QT1 
PRINT 51 
PRINT 39, Q 

CALL XTD (PR, PI, NEXP, KEXP) 

FACTOR = DY1-DY2 
DO 4 I = 1, N 
DO 4 J = 1, N 
PH = - (I-H+J-H)*CONST 
COSINE = COS(PH) 

SINE = SIN(PH) 

RE = PR(I, J)*C OSINE-PI(I, J)*SINE 
XM- PR(I, J)*SINE.+PI(I, J) "•'COSINE 
PR(I, J) = REDACTOR 
PI(I, J) = XM*FACTOR 

4 CONTINUE 

CALL SECOND (QT1) 

Q = QT1-QT 
QT = QT1 
PRINT 40, Q 
C 

C THIS COMPLETES THE CALCULATION OF THE UNBLOCKED 

C INTEGRAL. 

C 

c 

C CALCULATE THE CENTRAL BLOCKAGE CORRECTION. 
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C 

CALL, PHASE (0. , 0.,Z, NX, PH) 

TMAG = 2. 0*PIE*MAG(0. , 0. )*RB*RB 
RE = TMAG*COS(PH) 

XM= TMAG*SIN(PH) 

FAC1 = RB/DYl 
FAC 2 = RB/DY2 
DO 5 I = 1, N 
DO 5 J = 1, N 

F = 2.0*PIE/K*SQRT((FACl*(I-H)**2-f-(FAC2*(J-H))**2) 

F = F*F 

FACTOR = {(BES4*F+BES3)*F+BES2)*F+BES1 
PR(I, J) + PR{I, J)-RE*FACTOR 
PI(I, J) = PI(X, J)-XM*F ACTOR 

5 CONTINUE 

CALL SECOND (QTl) 

Q = QT1-QT 
QT = QTl 
PRINT 41, Q 
C 

C THIS COMPLETES THE CENTRAL BLOCKAGE CALCULATION. 

C 

C 

C CALCULATE THE SPAR CORRECTIONS. 

C JJ=1 FOR FIRST SPAR, AND 

C JJ=2 FOR SECOND SPAR. 

C 

DO 9 JJ = 1, 2 
DUP = DU1 

IF {DUP.GT. DU2) DUP = DU2 
FAC1 = COS{SIG(JJ)) 

FAC 2 = SIN(SIG(JJ)) 

DY = 2. *PIE/K/WAVENR/DUP 

CONST = PIE 

HK = K/2+1 

DO 6 I = 1, K 

X1P = (I-HK)*DY 

XI = X1P*FAC1 

X2 = X1P*FAC2 

CALL PHASE {XI, X2, Z, NX, PH) 

PH = PH-{I-HK)*CONST 
A{I) = MAG (XI, X2)#C OS(PH) 

B(I) = MAG (XI, X2)*SIN(PH) 

6 CONTINUE 

CALL FETR (A, B, K, 

C 



o o 
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APPLY PHASE TERM 

DO 7 X = 1, K 
PH = -COST*(I~HK) 

COSINE = COS(PH) 

SINE = SIN(PH) 

RE = A{I)*C OSINE-B(I)*SINE 
XM = A(I)*SINE-f-B(I)*COSINE 
A(I) = RE*W*DY 
B(I) * XM*W*DY 

7 CONTINUE 
C 

C THIS COMPLETES THE SPAR CORRECTION 

C IN THE PRIMED C ORDINATES. ' 

C 

C 

C APPLY SPAR CORRECTION.- 

C 

DO 8 I = 1, N 
DO 8 J = 1, N 
COSINE = COS(SIG(JJ)) 

SINE = SIN(SIG(JJ)) 

S = I-H 
T = J-H 

XV = S*DY/DYl*COSINE-T*DY/DY2*SINE+HK 
IV = XV 
IVP1 = IV+1 

RE = (A(IVP1)-A(IV))*XV+A(IVP1)-(A(IVP1)-A(IV))*IVP1 
XM = {B(IVP1)-B(IV))*XV+B(IVP1)-(B{IVP1}-B{IV))*IVP1 
PR(I, J) = PR(I, J)-RE 
PI(I, J) = PI(I, J)-XM 

8 CONTINUE 

9 CONTINUE 

CALL SECOND (QT1) 

Q = QT1-QT 
PRINT 42, Q 
C 

C THIS COMPLETES THE SPAR CORRECTION. 

C 

C 

C CALCULATE THE FAR FIELD GAIN. 

AGEO = P IE *DIAM *DIAM / 4 . 

FACTOR = WAVENR * WA VENR/PIE 
DO 10 I = 1, N 



152 


DO 10 J = 1, N 

PR(I, J) = FACTOR*(PR(I, J)**Z+PI(I, J)**2) 

10 CONTINUE 
GAIN = PR(H, H) 

F = 10. *ALOG10(GAIN) 

AF = LAMDA*LAMDA*GAIN/4. /PIE 
EFF = AF*100. /AGEO 
C 

C PRINT AXIAL GAIN PARAMETERS. 

C 

PRINT 29 
PRINT' 43, GAIN 
PRINT 44, AE 
PRINT 45, F 
PRINT 46, EFF 
DO 12 I = 1, N 
DO 12 J = 1, N 

IF (PR(I, J) .NE. 0.0) GO TO 11 
PR(I, J) = -40. 

GO TO 12 

11 PR(I, J) = 10. *ALOG10(PR(I, J))-F 

12 CONTINUE 
DO 13 I = 1, N 
DO 13 J = 1, N 

IF (PR(I, J) . LE. -40. ) PR(I, J) = -40. 

13 CONTINUE 
C 

C THIS COMPLETES THE CALCULATION OF 

C THE FAR FIELD .GAIN. 

C 

C 

C IN THE FOLLOWING A CONTAINS THE U1 PATTERN 

C AND B CONTAINS THE U2 PATTERN. 

C 

DO 14 I = 1, N 
A(I) = PR(I, H) 

B(I) = PR(H, I) 

D(I) = (I-H)*DU1*SCALE 
E(I) = (I-H)*DU2*SCALE 

14 CONTINUE 
C 

C CALCULATE BEAM WIDTHS. 

C 


D1 = DD1*1000 
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CALL WIDTH (A, 
CALL WIDTH (A, 
CALL WIDTH (A, 
CALL WIDTH (A, 
CALL WIDTH (A, 
CALL WIDTH (A, 
DT = DD2*1000. 
CALL WIDTH (B, 
CALL WIDTH (B, 
CALL WIDTH (B, 
CALL WIDTH (B, 
CALL WIDTH (B, 
CALL WIDTH (B, 


-03., 

Dl, 

D3, N) 

-10., 

Dl, 

D10, 

N) 

-15., 

Dl, 

D15, 

N) 

-20., 

Dl, 

D20, 

N) 

-25., 

Dl, 

D25, 

N) 

-30. , 

Dl, 

D30, 

N) 

-03., 

DT, 

P3, 

N) 

-10., 

DT, 

P10, 

N) 

-15., 

DT, 

P15, 

N> 

-20., 

DT, 

P20, 

N) 

-25., 

DT, 

P25, 

N) 

-30. , 

DT, 

P30, 

N) 


PRINT BEAM WIDTHS. 

PRINT 51 
PRINT 47 
PRINT 48 

PRINT 49, D3, DIO, D15, D20, DZ5, D30 
PRINT 48 

PRINT 50, P3, P10, P15, P20, P25, P30 

PRINT PATTERNS. 

PRINT 51 
PRINT 52 

PRINT 54, (A(I), I = 1, N) 

PRINT 51 
PRINT 53 

PRINT 54, (B(I), I = 1, N) 

PRINT PATTERN MAP. 

PRINT 29 
PRINT 55 
JJ = -15 

15 JJ = JJ-H6 
KK = JJ+15 
DO 16 J = 1, N 
II = N-J+l 

PRINT 54, (PR(I, II), I = JJ, KK) 

16 CONTINUE 
PRINT 29 

IF (KK .NE. N) GO TO 15 
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PLOT PATTERNS. 

XPLUS = H*DU1*SCALE 
YPLUS = H*DU2*SCALE 
XMINUS = -XPLUS 
YMINUS = -YPLUS 

CALL AXES (6. , XMINUS, XPLUS, -1., 4., -40., 0., 5., 5., 

132, 50) 

CALL PLOTITL (JX, 26, 4, 2, 0., -0.5) 

CALL PLOTITL (10, 3, 0, 2, -.395, 3.94) 

CALL PLOTITL (15, 3, 0, 2, -.395, 3.44) 

CALL PLOTITL (110, 3, 0, 2, -.395, 2.94) 

CALL PLOTITL (115, 3, 0, 2, -.395, 2.44) 

CALL PLOTITL (120, 3, 0, 2, -.395, 1.94) 

CALL PLOTITL (125, 3, 0, 2, -.395, 1.44) 

CALL PLOTITL (130, 3, 0, 2, -.395, 0.94) 

CALL PLOTITL (135, 3, 0, 2, -.395, 0.44) 

CALL PLOTITL (IDEC, 6, 0, 2, 4.625, 3.69) 

CALL PLOTITL (IPAT, 7, 0, 2, 4.805, 3.44) 

CALL PLOT (D, A, N, -8) 

CALL AXTERM (1) 

CALL AXES (6., YMINUS, YPLUS, -1., 4., -40., 0., 5., 5., 

132, 50) 

CALL PLOTITL (JX, 26, 4, 2, 0., -0.5) 

CALL PLOTITL (135, 3, 0, 2, -.395, 0.44) 

CALL PLOTITL (130, 3, 0, 2, -.395, 0.94) 

CALL PLOTITL (125, 3, 0, 2, -.395, 1.44) 

CALL PLOTITL (120, 3, 0, 2, -.395, 1.94) 

CALL PLOTITL (115, 3, 0, 2, -.395, 2.44) 

CALL PLOTITL (110, 3, 0, 2, -.395, 2.94) 

CALL PLOTITL (15, 3, 0, 2, -.395, 3.44) 

CALL PLOTITL (10, 3, 0, 2, -.395, 3.94) 

CALL PLOTITL (IPAT, 7, 0, 2, 4.805, 3.44) 

CALL PLOTITL (IPOL, 5, 0, 2, 4.95, 3.69) 

CALL PLOT {E* B, N, -8) 

CALL AXTERM (1) 

IF (YPLUS .GT, XPLUS) XPLUS = YPLUS 
XMINUS = -XPLUS 

PLOT MAP. 

CALL AXES (6., XMINUS, XPLUS, -1., 6., XMINUS, XPLUS, 5., 
15., 32, 1000) 

CALL PLOTITL (JX, 26, 4, 2, 0., -0.5) 

CALL PLOTITL (MM, 11, 0, 2, 4.5, 5.75) 
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CALL PLOTITL (I3DB, 4, 0, 2, 4. 86, 5. 50) 

CALL PLOTITL (110 DB, 4, 0, 2, 4. 86, 5.25) 

CALL PLOTITL (120 DB, 4, 0, 2, 4.86, 5.0) 

CALL PLOTITL (I30DB, 4, 0, 2, 4.86, 4.75) 

DO 27 JJ = 1, 4 

GO TO (17, 18, 19,20), JJ 

17 C = -3. 

GO TO 21 

18 C = -10. 

GO TO 21 

19 C = -20. 

GO TO 21 

20 C = -30. 

21 NN = 0 

KUP = N-l 

DO 25 I = 1, KUP 
XI = I 

DO 25 J = 1, KUP 
XJ = J 

IF ((FR(I, J) .GT. C) .AND. (PR(I, J+I) . LT. C) .OR. (PR(I, J) . LT. 
1C) .AND. (PR(I, J+1)'.GT. C)) GO TO 22 
GO TO 23 

22 XJ = J 
XJ1 = J+l 

XM = PR(I, J+l )-PR(I, J) 

P = PR(I, J)-XM*XJ 
NN = NN+1 

IF (NN . GE. 2045) GO TO 26 
YY(NN) = ((C-B)/XM~H)*DU2*SCALE 
XX (NN) = (XI- H)^DU1 -SCALE 

23 CONTINUE 

IF((PR(I, J) .GT. C) .AND. (PR(I+1, J) . LT. C) .OR. (PR(I, J) . LT. 
1C) .AND. (PR(I+1, J) .GT. C)) GO TO 24 
GO TO 25 

24 XII = XI+1 

XM = PR(I+1, J)-PR(I, J) 

B = PR(I, J)~XM*XI 
NN = NN+1 

IF (NN . GE. 2045) GO TO 26 
XX(NN) = ((C-B)/XM-H)*DU1*SCALE 
YY(NN) = (XJ-H)*DU2*SCALE 

25 CONTINUE 

26 CONTINUE 

CALL PLOT (XX, YY, NN, JJ) 

27 CONTINUE 
CALL AXTERM (1) 
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CALL AXTERM (0) 

C 

28 FORMAT (8F10. 3) 

29 FORMAT (1H1) 

30 FORMAT (//10X. *THE DIAMETER OF THE ANTENNA=*F6. 5* 

1FEET*/ / 1 0X*LAMDA=*F1 0. 3 * METERS*) 

31 FORMAT (//10X, *NEXP=*I3, 5X, *KEXP=*I3) 

32 FORMAT (//10X, *DELTA Ul=*F8. 6* DEGREES*5X*DELTA U2= 

1*F8. 6* DEGREES*) 

33 FORMAT (/ / 10X, *DELTA Ul=*Ei0.3* RADIANS*5X*DELTA U2 = 

1*E10. 3* RADIANS*) 

34 FORMAT (// 10X, *THE RADIUS OF THE BLOCKED AREA=*F8. 3* 
1INCHES*) 

35 FORMAT (//10X*THE ANGLE OF SPAR 1=*F6. 1* DEGREES*// 

110X*THE ANGLE OF SPAR 2=*F6. 1* DEGREES*) 

36 FORMAT (//10X*THE WIDTH OF THE SPARS=*F8. 6* INCHES*) 

37 FORMAT (//10X*THE EXTENT OF THE NONZERO APERTURE IN 
IT HE Y1 DIRECTION=*F6. 1* POINTS*) 

38 FORMAT (//10X*THE EXTENT OF THE NONZERO APERTURE IN 
1THE Y2 DIRECTION=*F6. 1* POINTS*) 

39 FORMAT (10X*THE TIME FOR INITILIZING ARRAY=*F16. 7* SECONDS*) 

40 FORMAT (10X*THE TIME FOR FFT=*F16.7* SECONDS*) 

41 FORMAT (10X*THE TIME FOR THE BLOCKAGE CORRECTION 
1*E16, 7* SECONDS*) 

42 FORMAT (10X*THE TIME FOR THE SPAR CORRECTION=*El6. 7* 

1 SECONDS*) 

43 FORMAT (/ / 10X*GAIN=*E16. 7) 

44 FORMAT (//10X*THE EFFECTIVE AREA=*F8. 3 * SQ. METERS*) 

45 FORMAT (/ / 10X*GAIN=*F8. 3* DB*) 

46 FORMAT (/ / 10X*EFFICIENCY=*F8. 3* PERCENT*) 

47 FORMAT (*WIDTHS (M-DEG) 3DB 10DB 15DB 20DB 25DB 30DB*) 

48 FORMAT (/) 

49 FORMAT (* U1 *6F7. 1) 

50 FORMAT (* U2 *6F7. 1) 

51 FORMAT {/////) 

52 FORMAT (* THE U1PATTERN FOLLOWS*) 

53 FORMAT (* THE U2PATTERN FOLLOWS*) 

54 FORMAT (2X1 6F7. 1) 

55 FORMAT (* PRINT MAP OF PATTERN*) 

END 
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SUBROUTINE XTD (PR, PI, NEXP, KEXP) 

SUBROUTINE XTD (PR, PI, NEXP, KEXP) 
DIMENSION PR(128, 128), PI(128, 128) 

DIMENSION A( 1024), B(1024) 

IF (NEXP . GE. KEXP) STOP 

N = 2**NEXP 

K = 2**KEXP 

JH = N+l 

DO 3 I = 1, N 

DO 1 J= 1, N 

A(J) = PR(J, I) 

1 B( J) = PI(J, I) 

DO 2 J = JH, K 
A(J) = 0. 0 

2 B(J) = 0. 0 

CALL FFTR (A, B, K, 0) 

DO 3 J - 1, N 
PR(J, I) = A(J) 

3 PI(J, I) = B(J) 

DO 6 I = 1, N 
DO 4 J = 1, N 
A(J) = PR(I, J) 

4 B(J) = PI(I, J) 

DO 5 J = JH, K 
A( J) = 0. 0 

5 B(J) = 0.0 

CALL FFTR (A, B, K, 0) 

DO 6 J = 1, N 
PR(I, J) = A(J) 

6 PI(I, J) = B(J) 

RETURN 

END 
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SUBROUTINE FFTR (X, Y, N, IFWD) 

SUBROUTINE FFTR (X, Y, N> IFWD) 

DIMENSION X(l), Y(l) 

COMMON/SCRFFT / F(1024), . G(1024)/FFF/S(5l2), C(512), IBR(1024) 
MASK = 00007777777777777777B 
NREM = NREM. AND. MASK 
IF (NREM-N) 1, 7, 1 

1 PI = 3. 1415926535897932384 
N2 = N/2 

RN = N$PIR = 2. *Pl/RN 
NEXP = ALOG10(RN)/ALOG10(2. ) 

NEXP1 = NEXP - 1 

C GENERATE SIN COS ARRAY S, C 

N4 = N/4$S(N4-M) = C(l) = l.$C{N4+l) = S(l) = 0.0 
N42 = N4+2$N22 = N2+2 
DO 2 I = 2, N4 
CC = COS(PIR*(I-l)) 

C(I) = CC$C(N22-I) = -CC$S(N42- 1) = CC 

2 S(N4+I) = CC 

C GENERATE BIT REVERSAL ARRAY IBR 
N1 = N- 1$N2 1 = N2-7 
DO 4 I = 1, N21, 8 

J = I-1$H = 0$ISH = -NEXP- 1$MASK = 1 
DO 3 L = 1, NEXP 
ISH = ISH+2 

NEW = LSHIFT(J, ISH) 

NEW = NEW. AND. MASK 
II = II-f-NEW 

3 MASK = LSHIFT (MASK, 1) 

ibr(i) = i+n 

4 IBR(I+N2) = IBR(I)+1 

IBR(3 ) = N4+1$IBR(5) = N4/2+l$IBR{7) = 3^4/241 

N7 = N-7 

DO 5 I = 9, N7, 8 

IDIF = IBR(I)-IBR(I-8) 

DO 5 J = 2, 6, 2 

5 IBR{I+J) = IBR(I4J-8)4IDIF 
DO 6 I = 1, Nl, 2 

JJ = N-I+l 

6 IBR(JJ) = N-IBR(I)+1 
C ADD UP 1ST ROW 

7 NREM = N 

IF (IFWD) 8j 8, 10 



159 


8 DO 9 1=1, N2 

IN2 = I+N2$F(I) = X(I)+X(IN2) 

G(I) = Y(I)+Y(IN2)$G(IN2) = Y{I)-Y(IN2) 

9 F(IN2) = X(I)-X(IN2) 

GO TO 12 

10 DO 11 I = 1, N2 

IN2 = I+N2$F{I) = X(I)+X{IN2) 

G(I) = - Y(X) - Y (IN2 )$G(IN2 ) = - Y(I)+Y(IN2) 

11 F(IN2) = X(I)-X(IN2) 

COMPUTE TREE 

12 JN = 2 

IF (N-4) 16, 16, 13 

13 DO 15 K - 2, NEXPl 

IE = N2/JN$NJN = N/JN$KK = -1 
DO 14 J = 1, JN 

KK = KK+2$1A = IBR(KK)$IS1 = (J-1)*NJN 
DO 14 I = 1, IE 
IS = IS1+I 
ISS = IS+IE 

SAVE = F(IS)$SAVER = G(IS) 

SAVEM = F(ISS)*C(IA)»G{ISS)*S(IA)$SAVES = F(XSS)*S(IAHG(ISS)*C(IA1) 
F(IS) = SAVEM+SAVE$G(IS) = SAVES+SAVER 
F(ISS) = - SAVEM+SAVE$G(ISS) = -SAVES+SAVER 

14 CONTINUE 

15 JN = JN*2 

C BIT REVERSE AND LAST ROW 

16 IF {IFWD) 17, 17, 19 

17 DO 18 I = 1, Nl, 2 

II = 1+1$ J = IBR(I)$JJ = IBR(II) 

SAVE = F(II)*C(J)-G(II)*S(J)$SAVER = F(II)*S(J)+G{II)*C{J) 

X(J) = F(I)+SAVE$X( JJ) = F (I) -SAVE 
Y( J) = G(I)+SAVER 

18 Y(JJ) = G(I)-SAVER 
RETURN 

19 DO 20 I = 1, Nl, 2 

II = 1+1$ J = IBR(I)$JJ = IBR(II) 

SAVE = F(II)*C(J)-G{II)*S{J)$SAVER = F(II)*S(J)+G(II)*C(J) 

X(J) = (F(I)+SAVE)/RN$X(JJ) = (F(I)-SAVE)/RN 
Y(J) = (-G(I)-SAVER)/RN 

20 Y(JJ) = ( - G(I)+SA VER) / RN 
RETURN 

END 
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SUBROUTINE WIDTH (A, C, DT, W, N) 

SUBROUTINE WIDTH (A, C, DT, W, N) 

C 

C THIS ROUTINE CALCULATES THE WIDTH AT 

C LEVEL C OF A WITH DISTANCE DT BETWEEN POINTS. 

C 

DIMENSION A(1Z 8) 

INTEGER H 
H = N/2+1 
DO 2 I = H, N 

IF ((A(I) .GT. C) .AND. (A(I+1) . LE '. C)) GOTO 1 
GO TO 2 

1 XM = A(I+1)~A(I)- 
B = A(I)-XM*I 

HW = ((C-B)/XM-H)*DT 
GO TO 3 

2 CONTINUE 

3 CONTINUE 
DO 5 I = 1, H 

IF ((A(I) . LE. C). AND. (A<I+I) .GT. C)) GO TO 4 
GO TO 5 

4 XM = A(I+1)-A(I) 

B = A(I)-XM*I 

HW2 = (H-(C-B)/XM)*DT 

5 CONTINUE 
W = HW+HW2 
RETURN 
END 
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FUNCTION MAG (X, Y) 

FUNCTION MAG (X, Y) 

C 

C THIS ROUTINE CALCULATES THE TAPERING 

C FUNCTION FOR A 2RO METER ANTENNA WITH A FOCAL 

C LENGTH FO. 

C GAIN IS THE GAIN OF THE FEED SYSTEM. 

C 

COMPLEX G 
REAL MAG 
G = (0. , 0. ) 

PIE = 3. 1415926535 
FO = 8. *12. *2. 54/100. 

RO = FO 

RSQ = X**Z+Y**2 
R = SORT (RSQ) 

IF (R.GT.RO) GO TO 1 
Z1 = RSQ/4. /FO 
PS = ATAN(R/(F0-Z1) 

SC = -3. 1415926535/2. 

IF (Y .GE. 0.0) XC = 3. 1415926535/2. 

IF (X .NE. 0. 0) XC = ATAN(Y/X) 

A = 1. 0 
B = 0. 66 

CALL FEED (GAIN, PS, XC, A, B, G) 

MAG = GAIN/4. /PIE/(FO+R*R/4. /FO)**2 
MAG = SORT (MAG) 

RETURN 
CONTINUE 
MAG = 0.0 
RETURN 
END 
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SUBROUTINE FEED (GAIN, TH, PH, A, B, G) 

SUBROUTINE FEED (GAIN, TH, PH, A, B, G) 

C 

C THIS ROUTINE CALCULATES THE DIRECTIVITY OF A 

C HORN ANTENNA IN DIRECTION TH=THETA 

C AND PH=PHI. THE DIMENSIONS OF THE HORN ARE 

C A/LAMDA AND B/LAMDA. 

C A 

COMPLEX CT1, CT2, CT3, CT4, G 

PI = 3. 1415926535 

G2 = 1.0-CABS(G)*CABS(G) 

X = PI*PI*Pl/2. /G 2 
U = SQRT(1. 25/A/A) 

S = SIN(TH) 

P COS(PH) 

T1 = PI*A*S*P 

FI = COS(Tl)/(Tl*Tl-PI*PI*. 25) 

FI = F1*F 1 
SP = SIN(PH) 

T2 = FI*B*S*SP 
F2 = 1. 0 

IF (T2 .NE. 0.0) F2 = SIN(T2)/T2 
F2 = F2*F2 
CT = COS(TH) 

Tl = 1. +U*CT 
CT 1 = T 1 *( 1 . 0, 0. 0) 

CT2 = CT1+G*CT1 
T3 = CABS(CT2) 

T4 = SP *SP #T 3 *T 3 

Tl = CT+U 

CT1 = Tl*(l. 0, 0. 0) 

T2 = CT-U 

CT2 = T2*(l. 0, 0. 0) 

CT3 = CTl+G*CT2 
T3 = CABS(CT3) 

F3 = T4+P*P*T3*T3 

GAIN = X *A*B *F 1 *F2 *F 3 / U 

RETURN 

END 



SUBROUTINE PHASE (XI, XZ, Z, NX, PH) 


SUBROUTINE PHASE (XI, X2, Z, NX, PH) 

PH=0. 0 

RETURN 

END 



APPENDIX B 


A. Introduction 

This Appendix gives the proofs of the astigmatic phase error 
properties. The phase error in the aperture plane under this model 
is given by Equation (25), 


r 2 2 2 2~| 2 

L(x) = j^or x^ + p x 1 x 2 " a x 2 + FK(xj + x^J/Rq . 

All of these properties take advantage of the fact that L(x) is even, i. e. 


L(x) = L(-x). 

B. Property 1 

If the illumination function E(x) is symmetric about the origin, i.e. 
E(x) = E(-x), then any contour map is symmetric about the origin for all 
feed positions and reflector errors. That is 


G(u) = G{-u). for all a, p, and T. 
Proof of Property 1 

The gain is given by Equation (7) as 


G(ur u 2 , a, b, r) - 4ttz p 


o m 


[ tti \ j k L ( x ) J k u ' x , | 2 
E(x) e — e — — dx| . 


Now by the assumptions of the property 


E(x) = E(-x) 

and 

L(x) = L(-x). 
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Making use of these assumptions and changing variables to y = - x, 
one has 


71 k 


G(u r V P ’ T) = 4ttz p. 


m 


| J E( Z ) e j k L(y > e j k <-2> ' X d Z | 


= Gf-u^ -u , a, p, r). 

This completes the proof of property 1 . 

G. Property 2 

If the illumination function E(x) is only a function- of r, the radius 
from the center of the aperture, when a contour map of the pattern is 
symmetric about the lines 


tan 2cp = — £■ 

2 o' 

for all feed positions. 

Proof of property 2 

In order to prove this property it is sufficient to show that 

g^ 7 , u z' ’ 6 ’ = u 2 *’ <*» 3. n 

and 

G^', , a, p,' D = Gfu^, -u ^ 7 , or, p, T) 

in a coordinate system u. 7 which is at an angle cp to the u coordinates. 
Only the first of these equalities will be shown. 

The u 7 coordinates are related to the u coordinates by the well 


known rotation 
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u ' = Ru, 


where 

r 

cos cp - sin cp 
sin cp cos cp 
Also let us define the x* coordinates by the same rotation 


R = 


Substituting u 7 into Equation (7), one has 


« 2 ', or, B,r) = 


T1 k 


-1 ,-T 


4ttz P. ' J 
o in 

x 


| f EM • e j k e j k[ S - 3 ‘ - <*S| 2 


Changing variables to x / in the above, one has 


H k 

G( u i » u 2 ’ a ’ ^ = 4ttz P7 


I f EMe jk{ “' (X l' 2 - X 2' 2) 


o in / 

x 


+r(x >Z + x /2 )]/R jk(u 'x '+u 'x '] 

1 i O iX 6 “ . / , / i ^ 

■ e dx x dx^ ] , 

where a' = a sec 2cp. 

Changing variables once again so that y = -Xp tt follows that: 


G{u 1 . u 2 , or, p; r) " 4ttz | j 


ue 1_ . r j em f t(y2 ' Xz,2) + r(y2 + X2 ' 2)5/r ° 


o m 


x 2 y 


jk{(-u 1 ')y + u 2 'x 2 } 


dy dx^l = G(-Uj , u^', a, g, f) 


This completes the proof of property 2, 
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D. Property 3 

If the illumination function E(x) is symmetric about the origin, i. e. 
E(x) = E(-x), then the gain is an even function of of, p, and F taken together. 
That is 


G{u r u . of, P, T) = Gfu^ u , -a, -p, -T). 


Proof of Property 3 

The gain function may be written by Equation (7) as: 


G(u r u 2 , or, 6, r) = | j E(x) e 


f 2 2 

jjcta X x + P x x x 2 " orx 2 


o in 


+ r( Xl 2 + * 2 2 )]/R 0 2 j k {»!*! + » 2 *2 5 2 

e dx] . 

The magnitude of a complex number is the magnitude of its complex con- 
jugant, so 


1 , k 2 jk «*“ x , 2 - » x , x 2 - (-"> x 2 2 

°<V V “• e ' r) P. I j E & e 


o m 


+ (-D ( Xj 2 + x/D/r/ jk C-UjXj -u 2 x 2 ] 2 

e dxj 

= G(-u lt -u^, - a , -p, -F) = Gfa^ u 2> -a, -p, ^T). 

The last equality follows from property 1. This completes the proof of 
property 3. 
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E. Property 4 

If the illumination function E(x) is only a function of r, then the 
angular coordinates which are aligned with the cp direction,, tan 2cp 

2of 

change roles when a and p change sign. That is 

G(Ui', u 2 7 , q;, p, D = G(u 2 /, ti^ 7 , -or, -P, F), 

where 


— — 


— 

— 


— — 

I / 






U 1 


cos cp 

- sin cp 


U 1 

/ 






U 2 


sin cp 

cos cp 


U 2 

— — 





^ 


Proof of Property 4 

From property 3, one has 


G{u^ ' , , a, 8, D = G{u J 7 , u^ 7 , -a, -R, -T). 

Then from property 5 to be proved subsequently, it follows that: 

G{Ui ' , u 2 7 , -o', -8, -D = g ( u 2 / . U 1 -a, -p, +F). 

This completes the proof of property 4. 

F. Property 5 

If the illumination function E(x) is only a function of r, then the 

angular coordinates which are aligned with the cp direction,, tan 2cp = -J^, 

2a 

change roles when T changes sign. That is 


Gfuj', , of, p, D = G(u 2 7 , u 1 / , or, p, -r). 


Proof of Property 5 

Referring to the proof of property 2, one has 
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T- k 2 n fria' (x / 2 - x/ 2 ) 

G(u /’ n z* ' ®* e ’ r) " 4 ttz P I j E(r)S 

o in z 


+r(xj ' 2 + x 2 ' 2 »/r o 2 ,, 1|2 

® d * 1 dx 2 l| • 

Once again observing that the magnitude of a complex number is the 
magnitude of its complex conjugant, one can rewrite the above equation 


H k 2 jk{o' / (x /2 - x /Z ) 

G( V- V- “• B ’ r) p l.f E(r)e 

o in , 


- r ( Xl ' 2 + X2 ' 2 ,)/R o 2 -jkfuj'xj' +V* 2 '1 , ,.Z 

e dx 1 dx^ j . 

Now changing variables to 


y l = " X 2 


y 2 = ~ X 1 ’ 


it follows that: 


G(u l ’ U 2 * V) = 4ttz P. i J E(r)e 


jk{o! / (y x 2 - y Z ) 


o m 


-r(yi 2 + y 2 2 »/R 0 2 jMa/y 2 + 

• d il 


= G(u 2 ' , Uj a, B, -f). 


This completes the proof of property 5. 



APPENDIX C 


This appendix gives two theorems which concern the frequency 

dependence of reflector antenna efficiency. Theorem 1 and 2 concern 

prime focus instruments with simple point source feed systems. Theorem 

1, however, is extended to apply to a Cassegrain antenna system. 

Theorem 1 requires as input data the design efficiency, T| , the 

measured efficiency at some frequency, T] 7 (f ) ± B , and an estimate 

- m m 

of the peak reflector phase error in radians, K(f }. From this infor- 

m 

mation Theorem 1 predicts both an upper and lower bound at another 

frequency, f . 

e 

Theorem 1 may be stated as follows: 

Theorem 1 

If the feed system of a prime focus reflector antenna at a frequency 

f is a scale model of the feed system at f , and 
e m 


Tl'(f )- 

m 


m )| 

m 1 


< B 


m 


and 


I 5(f , x) - 6{f . y)| < K(f ), all x and y, 

1 m — m — 1 — m — 

then the bounds 


f f 

We) 1 \ - p [\ - Tl ' <f m > j Gr) + \ 3 + p B m (r) • £ ° r a11 p 


m 


m 


and 
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f 2 f 2 

•«V> V P [V V<f m >j(r) forallP 


m 


m 


hold. The set S is defined as 


S = {z:z = U (~~'j - P (~'j U{x); 0 < x < K(f m ) }. 

m m 

The function U (x) is given by 


U(x) = cos (x) - 1. 

Proof of Theorem 1 

Recall from the text that the efficiency of a reflector antenna can be 
expanded into an infinite series in even powers of frequency. This series 
is given in Equation (31) as 


... V* „2n . 

W£) = L f A 2n 

n=0 

Denote by A 7 the estimate of A based upon the true efficiency at 

Lt 


f , i. e. 
m 


V =fz - V- 

m 


Let us define the estimate of the efficiency at frequency f as 


n (f ) = 11 + P A ' f 2 

e e o 2 e 


where P is a free parameter. Now one has 
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\ 7n z' p\ 

H(f ) - T1 (f ) = ; A f ~ T] - P[Tl(f ) - T1 1 Cf- ) . 
e e e L> 2n e o m o S 

m 


n=0 


Substituting of the series expression for Tj(f ) into the above equation and 
using straightforward manipulation, one obtains 


2 ® 


Z 2n 7 S & X 'T' 

A 2n f e + (1 - P > A 2 £ e - P (r) I A 2 n f 

m 


2n 
m ' 


n=2 


n=2 


Making use of a remainder function R[x] defined as R£x] = cos (x) + x /2 - 1, 
one has 


+f jI 

m 


t x 


E(x) E(x + t) |r[u]' -(1 - P) u 2 /2 - P J R j^*f”~ U jl d — d — ’ 

m e 


where 


u= 6{f , x) - 6 (f , x + t), 
e — e — — 

and 6 is the reflector phase error function. The proportionality factor, F, 
is defined in the text. 

The integral above is the mathematical error of approximation. It 
may be possible to actually perform the integral with 6 (x) being obtained 
from a STAIR program or mechanical measurement of the reflector surface. 
The result of this calculation would be a small correction to the efficiency 
estimate. The purpose here, however, is to bound this error. 
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Using the mean value theorem, one has 

t 2 

wy< v p c vvjdG 2 -) 

m 

and 

w e )>\-pc \-w m »Gr) + V nfs 

m 

where S is a set defined as 

f 2 f 

S = jz:z = R[u] - (1 - P) u 2 /2 - P J R u] ; 0 < u < K(f Q )}. 

me 

This definition for S is equivalent to the one given in the text. One must now 
take the measurement errors into account. Thus, 


K f e> S V P I V + \ SU ? S + P (r) 

m m 


B , for all P 
m 


and 


n(f)>11 - p Hu - *n # (f )] + Tt inf S - P B , for all P 

e “ o o m Vi y o Vf^m 

m m 

This completes the proof of Theorem 1. 

Theorem 2 is similar to Theorem 1. It uses a measurement at two 

frequencies to imply bounds at a third frequency. Theorem 2 may be stated 

as follows: 

Theorem 2 

If the feed system of a prime focus reflector antenna is scaled at 
frequencies f^, f^, and f^, and 
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j V^) - W L )| <b l 

- m 7 )\ 5 b ?- 


I 6{f , x) - 5{f , ^;)| < K(f ) all x and_£, 


and A ^ and A 17 are roots of the equations 


\ + V f i 2 + V l z - '"'‘V 


\ + V *** + V f 2 4 - ''««*>■ 


Then, 


Tj(f ) < H + PA " f 2 + QA // f 4 + IF.J B. + |F I B + H sup S, 
e — o 2 e 4e 1 1 1 1 * Z 1 2 o 


for all P and Q 


and 


T,(£ e ) > H o + PA/ f e 2 + QA/ £ e 4 - 1^1 Bj - |F 2 | B 2 + V nlS, 
for all P and Q, 

where 

f 'f 

S = {z:z = F 2 T[^ u j - Fj T[~ uj + T[u] - (1 - P) u 2 /2! 
e e 

+(1 - Q) u 4 /4!; 0 < u< K{f )} 

Proof- of Theorem 2 

The proof of Theorem 2 differs from that of Theorem 1 in that three 

terms of the cosine series are used instead of two. Define 1] (f ) as 

e e 
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T] (f ) = T| + P A " f 2 + Q A " f . 4 , 
e e o 2 e 4 e 


where P and Q are free parameters. Now one has 


- W * I f e 2n A 2n ' \ ' P V C ' Q V C 

n=0 

Let us define A* and A.' as the values A" and A." would have were there 
2 4 2 4 

no measurement errors. A sequence of straightforward manipulations yields 


T|(f e > - yy = P(A 2 - a 2 ') £ e 2 + Q(A 4 - a 4 ') f e 4 


+ (1 - P) A f 2 + (1 - Q) A f 4 + Y f 2n A • 
2 e 4 e L. e 2n 


n=3 


+ P(A,' - A") £ 2 + Q(A/ - A/) f 4 . . (Cl) 

2 2 e 4 4 e 

The first five terms are errors in approximation while the last two are 
caused by measurement errors. Clearly the summation term is 


CD 

) f 2n A = F f f f (x) f (x + t) T[ 6 (f , x) - 6 (f , x + t)] dx dT , 
i — i e 2n J J ” - — e — e — — — 

n=3 t x 

where 

T[x] = cos x - 1 + x 2 / 2! - x^/4!. 

One can now obtain expressions for A^ - A„ ' and A a - A ‘ . The definitions 

2 2 4 4 

for A^ 1 and A ' may be written as 


\ + A 2 ' £ l 2 + A 4' f l 4 = 1 « f l' 
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\ + A 2' f 2 2 + A 4' £ 2 4 = W 2 >- 


Rewriting this set of equations in matrix form, one has 


- 




— — 




CO 



V 


2 4 

£ i £ i 


A 2 


I £ : 2n A 2n 









n=3 




— 




+ 

CO 

2 4 

f f 
2 2 


V 


f 2 f 4 
2 2 


A 4 


, 2l1 A 

L f 2 A 2n 



M M 






_ n=3 


By inverting the square matrix above, it can be seen that and are 

A ' and A 7 respectively, plus error times. Thus 
c * * 


— 

***" 




CO 

V 

- A 2 




1 f l 2n A 2n 






n = 3 



1 




V 

- A 4 

< f lV 2 (f 2 2 - f l 2) 



ro 

^ f 2n A 
L 2 2n 






n=3 


and calculating at each error term individually one has 


V - A 2 = 


1 Tr 4 V 2n A . 4 r 2n -I 

f ,2 2 2 L f 2 L f l A 2n ' f l 1 f 2 A 2nJ 

' 12* (t 2 " 1 ' 


n=3 


n=3 


and 


~ A 4 = 


1 


(f f ) 2 (f 2 - f 1 
v 1 V ' 2 1 


[ f i 2 


V 

L 

n=3 


f 2n A - f 2 
2 A 2n f 2 


L 


, 2n i 

f l A 2nJ 


n=3 


Rewriting once again, one has 


A 2 - V 





F J J E{x) E(x + t) { T[6(f 2 »x) - 6(x + _t)] 


t x 
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f 

(~) TLMfjX) ~ 8(£ r x + T)]} dx dr 


(CZ) 


and 


A 4 ‘ A 4 , 2 , . 2 , 2 


f ? {£ ? ~ f f > 

Z Z 1 T X 


F j J E(x) E(x + t) |t[6(£ 2 ,x) - 6 + ^ 


- (j -J TH6{£ 1 ,3c +^r)Dj" cb£ djr- 


(C3) 


Now one needs an expression for the measurement errors A ' - A " 

c* c* 

and A 7 - A ^ . This may be accomplished by subtracting the defining 
relations for A ^ , and A^ from the defining equations for A ^ ' and A^ 7 , 

(A 2 ' - A 2 ") f : 2 + (A 4 r - A/) f x 4 = *n(f 1 ) - Tl / (f 1 ) = 


<V ' A 2 "> l Z + <V - A 4 > *Z = 1 '< £ 2> - 1 '' (f 2 ) = e 2’ 

where and are the measurement errors. Inverting this pair of 
equations one obtains 


1 


A ' - A " ={^£) - 

A 2 A 2 Vf^ f 2 f 2 


f l e o 

‘V f 2 - f 2 
2 *1 


(C4) 


and 


A ' _ A " = * 

A 4 A 4 „ 2 2 . 2 


'1 


V < £ 2 _ - V> £ 1 2 « f 2 2 - f l 2 > 


• (C5) 


i. e. 


Substitution of Equations (C2), (C3), (C4), and (C5) into Equation (Cl) gives 
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T,( £ e ) - H e <f e > - F j J EH E(x + T) { F 2 T [f u] - F, l[-± u] 

e e 

+ T[u] - (1 - P) u 2 / 2! + (1 - Q) u 4 /4!j dx dr + F x e - F 2 e 

where 


| Bj + | F 2 | B^, for all P and Q 

Tl(f ) > T| (f ) + 11 inf S - I F. I B. - I F,| B_, for all P and Q. 
e — ee o 1 1 1 1 '2' 2 

This completes the proof of Theorem 2. 

Theorem 1, Cassegrain Antenna Version 

The extension of Theorem 1 to the case of a Cassegrain antenna 

depends upon considering only the most significant terms in an expansion 

of the efficiency. A Cassegrain system in general has a slightly different 

illumination function at two different frequencies. The fractional change 


and 


u = 6(f ,x) -'8(f ,x + t), 

e — e — — 


F„ = 


, 2 , 2 ^ , 2 . 
f e {P f l " Q f e > 


2 f 2 (f 2 - f 2 ) 
2 '2 1 ’ 


C (P f Z 2 - ° O 

1 " S - 0 


Thus, since J e ^ J < B^ and | < B^ 

/ 

11(f e ) -\ (f e ) + T1 o SUpS + l F l 


and 
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in the illumination at the two frequencies will be given by e{x), i. e. 

E(f , x) = E(f , x + t) [1 + e(x)], 
e — m — — 

where f is the frequency of estimation and f is the frequency of 
e m 

measurement. * In addition the feed system of a Cassegrain antenna 
system usually introduces a small phase error. This error will be 
denoted by k-t(x) . 

With these definitions and the assumptions listed in the text the 
stated bounds can be derived. These added assumptions are as follows: 

1. Integrals involving the sine of the reflector phase 
error differences between two points in the aperture 
are negligible. The Ruze model for the reflector error 
predicts that the expectation of these integrals are 
identically zero. In addition the Ruze model will be 
used explicitly to evaluate one term. 

2. Fourth order terms and above in e(x), -t{x), and 6{x) 

4 

can be eliminated. The 6 (x) terms are retained, 
however. 

The efficiency at frequency f is 

T)(f ) = F f | E(f , x) E(f , x + T ) cos [S + T ] dx d T 

G J J G " G *“ *""* G G ““ 

T X 

where E{f^, x) is the magnitude of the electric field in the aperture plane 
at the frequency f . The functions S and T are 

G G G 
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S - Z(£ , x) - -t»{f , x + t) 

6 6 — 0 — — 

and 

T = 6 (f , x) - 6 (f , x + t). 
e e — e — — 

Expanding the cosine of a sum factor, the above expression 
becomes 

T|(f ) = F f [ E(f ,x) E(f , x + t) cos S cos T dx dx 
e J J e~ e — — e e~~ — 

_X X 

— F I r E{f , x) E(f , x + t) sin S sin T dx dx. 

<1 J e e — — e e — — 

.I 

The second integral in this expression is zero under assumption 2 above. 

It is felt that this term should be very small for any reasonable phase 
function 6 . 

One can now expand the cosine of T term in the first integral so 

0 

that 

Tl(f ) = H <f ) " F f f E(f , x) E(f , x + x) T 2 /2! dx dx 
e oe jJ e— e — — e — — 

JT X 

+ F f r E(f , x) E(f , X + x) R[T 3 dx dx 
JJ e — e — — e — 

x x 

where R[x] is defined as 

2 

R[x] = cos x - 1 + x / 2! . 

The cos S term has been eliminated from the second integral by assumption 
e 

2. The leading term is by definition the design efficiency, H (f ), the 

o e 
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efficiency without phase errors in the primary reflector. 

Inserting the free parameter P and tahing advantage of the fact 
that the reflector errors scale, one has 

.£ 2 

1\(f ) = n (f ) - p (-r- ) F \ I E(£ , x) E(£ , x + t) T i 12’. dx dr 
e oe Vf y J J e— e — — m — — 

m 

T x 


+ fJ j E(f , x) E{f e> x + t) {{P - 1) T g e /2! + R[T e ]j- dx dr. 


T X 


Substituting the fractional change in illumination into this expression, 
it follows that 


T|(f ) = 1) (f ) - ~P(~) F f f E{f ,x) E(f , x + T ) T 2 /2i dx d T 
' e o e Vf > J J m — m — — m — — 

m 


- T X 


-P F j J [e(x) + e{x + _r) + e(x) e(x + t)] ^(f^, x) Eff^, x + t) 

m 

T x 


T 2 / 2! dx dr 

m — — 

+F [ f E(f , x) E(f . X + t) {(P - 1) T 2 /,2! + R[T dx dr. 

J J S © V Q J 

T X 

The second integral is now evaluated by the Ruze model after fourth 
and higher order terms are dropped so that the expansion becomes 
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Tl(f ) = H (f ) - P > F f f E(f , x) E(f , x + t) T 2 / 2! dx dr 
e oe vf y JJ m — m ^ — m — — 

m 


T x 



m 


+ fJ j E(f e , x) E(f e , x + t) {(P - 1) T Z /2! + R£T ]j- dx dr. 

T X 

The task now is to determine an expression for the first integral 

above in terms of the efficiency at the measurement frequency. This can 

be done by taking the difference between the true efficiency, Tl (f ), at f 

m m 

and the design efficiency at f , T1 (f ): 

m o m 


T|(f ) - 1) (f ) = F I I E(f , x) E{f , x + t) cos S cos T 
m om J J m — m — — m m 

x 

- f[ f E(f , x) E(f , x + t) COS S . 

J J m— m — — rn 

_T X 

In the first integral the sin T terms were eliminated by assumption 1 

m 

as before. Expanding the cosine terms in a power series and applying 
assumption 2, one obtains 

*n'(f ) - II (f ) = - F f f E(f , x) E(f , X + T) T 2 /2! dx dr 

m om JJ m — m — — m — — 

T x 


+ F 



E(f , x) E(f , x + T ) R[T ] dx d T . 
m — m — — m — — 


t x 
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Solving for the first integral and substituting the result into Equation 
one has 

A. 2 

m e ) = V f e> - P (rJ C W - Wm'l 


+ P ° <0 - W ] 


i 

- F I J E t f m ' *> E <w ^ + 1> { p (r) R[T m ] ; * 

m 

T X 


+ F j j E(f e , x) E(f g , X + JT) |{P - 1) T e /2! + R[T e ] j dx dr. 

T_ X 

Once again making use of the fractional change in illumination and elimi- 
nating fourth and.higher terms, one has 

2 


■<i<y = yy - p vjv [ W ■ '"m 11 

m 


+ “ < f a> W - W] 


+ r j j E<f m . *> E(f m , * + t> (<p - 1) t// 2! + R[T e ] - P(f ) R[Tj} 


T X 


dx d_r. 

The bounds given in the text now follow from the mean value theorem: 


w e > < yy - p (r) cy £ m > - vtf m n + p (r) B 

m nr 
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+ a (f e ) - H 0 (f e )] + \(f m ) sup S, for all P 


and 


f 2 f 2 

^ ft N 


[ V £ m > - - p (r) e 

in 


m 


+ a < f e ) £ V f m> ~ \ (f e^ + V £ m > inf S > for a11 p - 


The set S is given by 

2 


f 

S = jz*.z = (P - 1) u^/2! + R[u] - P (-f~J R[u]; 0 £ u < K(f )j 

m m 


This definition for S is equivalent to the one given in the text. 
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